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Abstract 



The main purpose of this master thesis is to study the LS-sequences of points 
introduced by Carbone in [2] and find two generalizations of them to the unit 
square. Here we also present a new algorithm proposed by the sanie author 
in [3] and we implement it in order to have a graphical description of these 
sequences. 

Chapter 1 includes a collection of results concerning the uniform ditribution 
theory and the discrepancy (we refer to [1] and [8] for a complete survey on 
the matter). 

In Chapter 2 we focuse our attention on the LS'-sequences of partitions and 
of points in the unit interval, giving particular attention to the ordering of 
the points "à la van der Corput "and finding a way to compute them related 
to the digit expansion of naturai numbers in base L + S. 
In Chapter 3 we move on the unit square where we find two generalizations of 
the LS'-sequences following the historical development of the van der Corput 
sequence in the multidimensional case. This is the reason why we cali the 
first sequence LS- sequence of points à la van der Corput-Hammersley and 
the second one LS'-sequence à la Halton. 



Introduzione 



Da lungo tempo lo studio della Teoria dell' Uniforme Distribuzione suscita 
grande interesse da parte di molti matematici. 

Gli albori di questa teoria risalgono agli inizi del '900 e precisamente l'anno 
di nascita viene fatto coincidere con il 1916, quando H. Weyl pubblica il la- 
voro dal titolo Uber die Gleichverteilung von Zahlen mod. Eins [H], in cui 
generalizza i risultati ottenuti dal matematico Kronecker sull'approssimazio- 
ne di equazioni diofantee. 

Da quella data molti passi in avanti sono stati compiuti ed il campo di in- 
dagine si è esteso dalla Teoria dell'Approssimazione a rami della matematica 
che vanno dall'Analisi Funzionale alla Teoria dei Numeri, dalla Teoria della 
Probabilità all'Algebra Topologica. 

Lo scopo principale di questa tesi è stato quello di estendere al quadrato 
unitario le LS-successioni di punti con L, S G N soddisfacenti un'opportuna 
condizione, introdotte, insieme alle LS- successioni di partizioni, da I. Car- 
bone in [2] nell'intervallo [0, 1[. 

In particolare, questo lavoro si articola in tre capitoli così ripartiti. 

Nel primo capitolo si è scelto di presentare i concetti di base ed i criteri 
classici della Teoria dell'Uniforme Distribuzione modulo 1, necessari per la 
comprensione dei capitoli successivi e che permettono di stabilire quando una 
successione di numeri reali è uniformemente distribuita nell'intervallo unita- 
rio, seguendo l'impostazione di [8]. Tra questi, uno dei risultati più rilevanti, 
in quanto si tratta di un criterio sufficiente per determinare classi di succes- 
sioni uniformemente distribuite, è il Criterio di Weyl, di cui abbiamo preso 
anche in esame alcune tra le sue molteplici applicazioni. 
Queste nozioni sono state poi estese, in maniera naturale, a successioni a 
valori in W. 

Infine, si è scelto di presentare i risultati esposti in precedenza da un punto 
di vista quantitativo e abbiamo introdotto la quantità detta discrepanza che 
ci permette di distinguere tra "buone" e "cattive" successioni di punti uni- 
formemente distribuite. Infatti, tra le successioni uniformemente distribuite 
modulo 1 ve ne sono alcune che presentano una buona distribuzione, al con- 



trario di altre che sono soltanto uniformemente distribuite. 

Il Capitolo 2 è caratterizzato dall'estensione della teoria dell'uniforme distri- 
buzione alle successioni di partizioni, concetto introdotto dal matematico S. 
Kakutani in [7j nel 1976. In particolare, abbiamo preso in esame la costruzio- 
ne di successioni di partizioni con la tecnica degli «-raffinamenti introdotti 
dallo stesso Kakutani e più in generale il metodo dei p-raffinamenti introdotti 
da A. Volcic in [12] e abbiamo presentato la classe delle LS-successioni di 
partizioni in [0, 1[ che fanno uso di tale costruzione. 

Abbiamo visto che esiste uno stretto rapporto tra le successioni di partizioni 
e le successioni di punti ad esse associate e abbiamo presentato due metodi, 
introdotti da Carbone in [2] e [3], per la determinazione di queste ultime, 
cercando di mettere in luce come il secondo metodo sia più efficiente rispetto 
al primo. 

Abbiamo anche dato la definizione di discrepanza di successioni di partizioni, 
e presentato le stime per la discrepanza delle LS'-successioni di partizioni e 
anche per quelle di punti ad esse associate, dimostrate in [2]. 
In questo capitolo, inoltre, abbiamo dato grande spazio agli esempi in modo 
da rendere la trattazione degli argomenti più chiara e a questo scopo molti 
degli esempi sono stati corredati di grafici. Infine, abbiamo concluso questa 
seconda parte con la presentazione di un algoritmo che descrive la procedura 
"alla van der Corpuf'per determinare le LS'-successioni di punti. Quest'ul- 
timo ci ha permesso di implementare i risultati descritti precedentemente. 

Nel terzo capitolo abbiamo dapprima considerato gli esempi classici di suc- 
cessioni uniformemente distribuite sul quadrato di lato unitario, cioè le suc- 
cessioni di J. C. van der Corput [11], J. M. Hammersley [6] e J. H. Halton 
[5]. Trattandosi di successioni aventi discrepanza asintotica ottimale (dette 
a bassa discrepanza) abbiamo descritto brevemente il Metodo Quasi-Monte 
Carlo, che fa uso di successioni di questo tipo. 

Anche una sottofamiglia numerabile delle LS-successioni di punti (e di parti- 
zioni) è a bassa discrepanza. Pertanto, in maniera analoga a quanto fatto da 
van der Corput, nel passaggio dalla dimensione 1 alla dimensione 2, costruia- 
mo una prima generalizzazione delle LS-successioni di punti sul quadrato che 
abbiamo chiamato LS-successioni "alla van der Corput". I risultati ottenuti 
fanno ben sperare che si tratti di una successione uniformemente distribuita. 
In seguito, siamo passati alla trattazione vera e propria di successioni (infi- 
nite) di tipo LS sul quadrato, costruendo coppie di LS-successioni di punti 
che abbiamo chiamato LS-successioni "alla Halton". Questi ultimi studi han- 
no portato alla luce risultati molto variegati tra loro. Dai grafici di queste 
successioni, infatti, si possono osservare distribuzioni di punti abbastanza 
soddisfacenti ed altre che presentano un comportamento alquanto bizzarro. 



Tra queste ultime, quella che ha suscitato il maggiore interesse è stata la 
successione ottenuta considerando le coppie L — S — ledL — 4, S — 1, il 
cui grafico risulta essere il più irregolare come vedremo in dettaglio nel corso 
del capitolo. 

La tesi si conclude con la congettura che le LS'-successioni "alla van der Cor- 
put" sono uniformemente distribuite. Si tratta di una congettura, dal mo- 
mento che nessun risultato teorico avalla ancora la nostra tesi, ma i risultati 
numerici e grafici ottenuti in questa tesi rendono questo risultato plausibile. 



1 



Capitolo 1 



Successioni di punti 
uniformemente distribuite 
modulo 1 



Lo scopo di questo primo capitolo è quello di presentare i concetti di base 
della Teoria dell' Uniforme Distribuzione, necessari per la comprensione dei 
capitoli successivi. Inizieremo con alcune definizioni ed alcuni risultati classici 
sulle successioni di punti uniformemente distribuite, forniremo la definizione 
di discrepanza di una successione di punti e presenteremo alcuni esempi ben 
noti in letteratura, seguendo l'impostazione classica di [8]. 

1.1 Definizioni e nozioni preliminari 

Sia (x n ) n£ N una successione di numeri reali. Ogni x n può essere scomposto 
nella somma della sua parte intera [x n \ = sup{m G Z | m < x n } e della 
sua parte frazionaria {x n } = x n — \_x n \. Dato che la parte frazionaria di x n 
coincide proprio con il valore di x n modulo 1, lo studio di successioni (x n ) 
di numeri reali modulo 1 si riconduce allo studio delle successioni delle parti 
frazionarie di quest'ultime nell'intervallo [0, 1[. Se la successione è a valori 
nell'intervallo unitario, allora essa coincide con la successione delle sue parti 
frazionarie. 

Definizione 1.1.1. Una successione (x n ) ng N di numeri reali si dice unifor- 
memente distribuita modulo 1 (u.d. mod 1) se per ogni coppia di numeri reali 
a e b, con < a < ò < 1, si ha 





n=l 
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dove X[a,b[ denota la funzione caratteristica dell'intervallo [a, b[. 

Una definizione equivalente si ottiene se consideriamo intervalli del tipo 
[0, c[ , con < c < l. Infatti, vale il seguente risultato. 

Proposizione 1.1.2. Una successione (x n ) ne ^ di numeri reali è u.d. mod 1 
se e solo se 



1 N 



^5>[M(K}) = c (1.2) 

71=1 

per ogni numero reale c in ]0, 1]. 



Dimostrazione. La dimostrazione è immediata. Infatti, se (x n ) è una suc- 
cessione di numeri reali e se [a, b[ è un sottointervallo dell'intervallo unitario, 
allora 



1 N 



n=l 

1 N 1 N 

= n £ WKK}) - x E X[0,a[({^n}) = 6 - a • 

n=l n=l 

Dunque la successione è u.d. mod 1 in [a, b[. 

Il viceversa è ovvio. □ 

Per brevità di notazione, nel seguito porremo I = [0, 1[. 
Vediamo adesso un primo risultato che prende il nome dal matematico te- 
desco H. Weyl il quale, nel 1914, fu il primo a dimostrare questo importan- 
te criterio che fornisce una condizione necessaria e sufficiente affinché una 
successione sia u.d. mod 1. 

Teorema 1.1.3. Una successione (a; n )„ e N è u.d. mod 1 se e soltanto se per 
ogni funzione f definita su I, continua e a valori reali, vale 



1 N r 1 

iV->-0O iv L — ' In 
n=l JiJ 



Dimostrazione. Dimostriamo che la condizione è necessaria supponendo dap- 
prima che la / sia una funzione semplice definita sull'intervallo unitario. 
Sia dunque (x n ) u.d. mod 1 e la / definita da 

fc-i 

i=0 

dove = a < ai < . . . < a k = 1 e Ci E M. per ogni i — 0, . . . , k. 
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Da questa relazione e dalla (1.1) segue che 

N 



N k-1 



n=l 



AT-s>oo jV 



re=l i=0 



k— 1 fc— 1 »i »i 

^Ci(a m - a,) = / CiX[ aii(Ii+l[ (a;)da; = / 



f(x)dx . 



Sia ora / una funzione reale continua su /. Si tratta dunque di una funzione 
integrabile secondo Riemann e quindi, fissato e > 0, esistono due funzioni 
semplici f%, f 2 tali che 



ed inoltre 



Pertanto 



fi(x) < f{x) < h{x) VxG/ 

/ [h{x)-h{x)\dx<e. 
Jo 

fi(x)dx< ì f 2 (x)dx < / fi(x)dx + e . 



Da questa catena di disuguaglianze e dalla (1.1) si ha che 

f 1 f 1 1 N 

/ f(x)dx -e< fx(x)dx = lini — V fi({x n }) 
Jo Jo N ~>°° N ^ 

i N ii! i N 

AT-s-noiV — iV->ooiV — ' jV-S-oo iV z — ' 



f 2 (x)dx < / /i(x)da: + £< / /(a^Gfe + e. 



Dall'arbitrarietà di e segue la tesi. 

Viceversa, sia (x n )„ g N una successione di numeri reali che verifica la (1.3) 
per ogni funzione continua in I a valori reali. Sia [a, 6[c / un arbitrario 
sottointervallo di /. Fissato e > 0, esistono due funzioni continue g\ e g 2 tali 
che 

9i(x) < X[a,b[(x) < g 2 (x) \/xeI 

ed inoltre 

ri 

[g 2 (x) - gi{x)]dx < e . 
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Allora 

b — a — e 



X[a,b[( x ) dx - 8 ^ / 92(x)dx-e< / gi(x)dx 



1 N 1 N 



n=l 

iV 



n=l 



n=l 



7V->oo TV 



n=l 



g 2 (x)dx < / X[ a ,6[(a;)c?x + 



6 — a + e . 



Essendo e arbitrariamente piccolo, segue la (1.1) e quindi la tesi. 



□ 



Conseguenza del precedente teorema è il seguente 

Corollario 1.1.4. Una successione (x n ) ne ji è u.d. mod 1 se e solo se per 
ogni funzione f integrabile secondo Riemann su I vale la (1.3). 

Dimostrazione. Per la condizione sufficiente occorre solo osservare che ogni 
funzione continua è integrabile secondo Riemann e quindi applicare il teorema 
precedente. 

La parte non banale della dimostrazione consiste nel provare la condizione 
necessaria. Procediamo in maniera analoga a quanto fatto nella prima parte 
del Teorema 1.1.3 Supponiamo che / sia una funzione integrabile secondo 
Riemann e che valga la condizione (1.1). Allora, fissato e > 0, esistono due 
funzioni semplici hi, h 2 tali che 

hi{x) < f(x) < h 2 {x) Vx e I 



[h 2 {x) — hi(x)]dx < 



Allora esiste N tale che per ogni N > N si ha 



1 N f 1 
— ^/i fc ({x n }) - / h k (x)dx 

n=l ^° 



< - per k — 1, 2 



perchè la (1.3) vale per funzioni semplici. Di conseguenza 

e 1 N 1 N 

hi(x)dx-- < jj^hidxn}) < — ^2f{(x n }) 



n=l 
JV 



n=l 



i r 1 

< ^J2 h ^ x ^< / ^ 

n=l ^° 



[x)dx + 
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Dalla condizione 

■i ,i ,i 

h\{x)dx < / f(x)dx < / h 2 (x)dx 
Jo Jo 

segue la tesi. □ 



Vediamo ora una generalizzazione del Teorema 1.1.3 al caso di funzioni 
periodiche a valori complessi, che ci sarà utile nel paragrafo successivo. 

Corollario 1.1.5. Una successione (x n ) n ^ è u.d. mod 1 se e solo se per 
ogni funzione continua f : IR — > C periodica di periodo unitario si ha 

i N r 1 

lim ^ /(*„)= / f(x)dx. (1.4) 

N^oo TV z — ' /n 

n=l JU 

Dimostrazione. Supponiamo che / : R — > C sia continua. Allora, applicando 
il Teorema 1.1. 3| separatamente alla parte reale ed alla parte immaginaria 
di /, si ottiene che la (1.3) vale anche per funzioni a valori complessi. Per 
ottenere la ( |1.4 ) è sufficiente osservare che la periodicità della / implica che 
f({ x n}) = f(x n ). 

Per dimostrare l'altra implicazione basta richiedere che le funzioni g\ e g 2 , 
utilizzate nella dimostrazione del Teorema 1.1.3 soddisfino l'ulteriore con- 
dizione <7i(0) = gi(l) e (72(0) = (72(1), cosicché entrambe possono essere 
prolungate su tutto IR per periodicità e soddisfano la (1.4). □ 

La prossima proposizione rappresenta un metodo efficace per generare 
nuove successioni u.d. a partire da una successione di cui è nota la proprietà 
di uniforme distribuzione. 

Proposizione 1.1.6. Se (x n )neN è una successione u.d. mod 1 e se (y n )neH 
è una successione tale che lim^oo \x n — y n \ = 0, allora (y n ) n6 N è u.d. mod 
1. 

Dimostrazione. Fissiamo a G ]0, 1[ e e > tale che e < min{a, 1 — a}. Allora 
dall'ipotesi lim^oo \x n — y n \ = segue che, per n sufficientemente grande, 

se {x n } < a — e allora {y n } < a 

e se {y n } < a segue che {x n } < a + e . 



Pertanto 



l n 1 N 

a - e = lim — V X[o,a-e[{{%n}) < lim — V X[o, a [({z/ n }) 

n=l n=l 

1 N 1 N 

- } ìm M X[0,a[({yn}) < l r im 22 X[0,a+e[({Xn}) =a + € . 
N— s>ooTV ^ — ' N~>oo TV — ' 

n=l ra=l 
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Nuovamente dall'arbitrarietà di e segue la tesi. 



□ 



1.2 II Criterio di Weyl 

In questo paragrafo tratteremo l'importante criterio di Weyl. Esso prende 

, dove h è un intero non nul- 
lo. Dimostreremo che queste funzioni soddisfano le ipotesi del Corollario 



in esame le funzioni della forma f(x) 



1.1.5[ e pertanto esse forniscono un criterio sufficiente a determinare se una 
successione è u.d.. 

Teorema 1.2.1 (Criterio di Weyl). La successione (x ra ) ne N è u.d. mod 1 se 
e solo se 



lim 



1 N 

v 



Jl-KÌhx n 







1.5) 



71=1 



per ogni intero h ^ 0. 

Dimostrazione. Sia (x n ) una successione u.d. modi. Allora, per il Corollario 



1.1. 5| e per la periodicità della funzione integranda, si ha 

x=l 



lim . 

Ar->oo N 



1 N 
V 



^2-TTÌhXn 



71=1 



,,2-Kihx 



-,2iTÌhx 



dx 







2nih 



. 



x=0 



Viceversa, supponiamo che la successione verifichi la condizione (1.5) e che 
la / sia una generica funzione continua su [0,1], con /(0) = /(l). Per il 
teorema di approssimazione di Weierstrass, fissato e > esiste un polinomio 
trigonometrico ip(x), cioè una combinazione lineare di funzioni del tipo e 2mhx , 
tale che 

sup \f(x) - ip(x)\ < - . 
xe[o,i] ó 

Allora 



+ 



N 



< 


L 







{f{x) — ip(x))dx 



1 1 N x N 

^ 71=1 71 = 1 



+ 



< £ 



poiché, per la (1.5), il secondo addendo è maggiorato da | per iV sufficiente- 
mente grande. □ 
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Applicazioni del Criterio di Weyl 



Volgiamo ora la nostra attenzione ad alcune applicazioni del criterio di Weyl. 
Nella prima parte di questo paragrafo prendiamo in esame una successione 
numerica speciale. 

Esempio 1.2.2. Sia 9 un numero irrazionale. Allora la successione (n9) n ^ 
è u.d. mod 1. Infatti 



N / N 

1 \ ^ 2-KÌhnO _"_ [ \ ^ 2-KÌhnd 

N ^ N\^> 

n=l \n=0 



N \ 1 — e 2irihe 



1 

N 



l — e 2mh6(N+l) _ ]_ _|_ e 27r * 



he 



1 -e 



2irih8 



1 



2mh9\ 



1 

N 



N 

l _ e 2nihN9 



1-e 



2mhNe 



1-e 



2nih9 



l — e 2nih9 

A questo punto osserviamo che, posto a = irh9, si ha: 



Il -e 



2ia\ 



|1 — cos(2a) — % sin(2a) 



1 - cos(2a)) 2 + sin 2 (2a) 



a/2(1 — cos(2a)) = y 2(1 — cos 2 a + sin 2 a) = 2| sina| 



Allora vale la (1.5) perchè 



1 N 



,,2irihn9 



n=l 



1 

N 



n 2irihNe 



2 sm(irhe) 



< 



1 



iV|sin(7r^)| 



è una quantità infinitesima per N che tende ad infinito. Dal Teorema |1.2.1 
segue che la successione {nO) è u.d. mod 1. 

Il fatto che la successione (n8) ne ?q sia u.d. mod 1 se 9 è un numero irrazionale 
fu stabilito in maniera indipendente da P. Bohl in pQ, W. Sierpinski in [U] e 
[TU] e Weyl in [J3] nel 1909-1910. L'importanza di questo risultato risiede nel 
fatto che esso segna l'inizio della teoria dell'uniforme distribuzione. Infatti 
la successione (n9) fu il primo esempio di successione u.d. mod 1 fornito 
dallo stesso Weyl ed esso perfeziona l'enunciato di un teorema di Kronecker, 
il quale asserisce che la successione (nO) è densa in [0, 1]. 



Presentiamo ora un risultato che è una conseguenza del criterio di Weyl 
e che ci consentirà di determinare molte classi di successioni u.d.. 
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Teorema 1.2.3. Sia (f(n)) ne ^ una successione di numeri reali tale che 
la successione Af(n) = f(n + 1) — f(n), n = 1,2,..., sia monotona. 
Supponiamo inoltre che 



lim A/(ra) = e lim n\Af(n)\ = +oo . 
Allora la successione (f(n)) è u.d. mod 1. 



1.6) 



Dimostrazione. Fissiamo due numeri reali u e v. Allora, integrando per parti 
47T 2 Jq~ v (u — v — w)e 2mw dw ed effettuando le opportune maggiorazioni si ha 
che 

|e 2 ™ - e 2 ™ - 2m(u - v)e 2niv \ = \e 27riv \ \ e 2 ^ u ~ v ^ - 1 - 2ni(u - v)\ = 



= 4tt 2 



4ir z 



(u-v- w)e 2mw dw 



w 



2\ 1 



U — VÌW — 



<4tt 2 



4tt 2 



\u — v — wìdw = 



(u — v)' 



2ix 2 (u-vf 



w=0 



Ponendo nella disuguaglianza precedente u = hf(n + 1) e v — hf(n), con 
h E Z, h 7^ 0, si ottiene 

| e 2*ifc/(n+i) _ e 2^/(n) _ 27 rihAf(n)e w(n) | < 2vr 2 /i 2 (A/(n)) 2 (n > 1) . 
Dividendo ambo i membri della disuguaglianza per |A/(n)| si ha 



D 2nihf(n+1) 2mhf(n) 



Af(n) Af(n) 



Da questa relazione segue che 



27rihe 2mhf{n) 



< 2ix 2 h 2 \Af(n) 



2wihf(n+l) 2mhf(n) 



A/(n + l) A/(n) 



2nihe 2mhf{n) 



< 



,2irihf(n+l) 2irihf(n+l) 



A/(n + l) A/(n) 
1 

< 



+ 



2irihf(n+l) 2mhf{n) 



Af(n) Af(n) 
1 



- 2mhe 2 * ihf ( n) 



A/(n + l) A/(n) 



+ 27r 2 /i 2 |A/(n)| 
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per n > 1. 

Dunque otteniamo la stima 



N-l 



} 2irihf(n) 



n=l 



N-l 



[2me 2mh ^ - 



e 2irihf(n+l) e 2nihf{n)\ e 2itihf{N) e 2nihf(l) 



+ 



+ 



N-l 

SE 

n=l 

N-l 

- E 



2nie 27rihf{n) - 



A/(n+l) A/(n) ; A/(iV) A/(l) 

e 27TÌhf(l) 



n=l 



e 2irihf(n+l) g2wihf(n) 

A/(n + l) + "A/H" 

AT-l 





e 2irihf{N) 




+ 


Af(N) 


+ 



A/(n+l) A/(n) 



A/(l) 



n=l 



Da ciò, sfruttando l'ipotesi di monotonia si ottiene 



N-l 

E 



n=l 



< 



1 



1 



1 



1 



N-l 



+7r|fc| £ |A/( 



2n\h\ y\Af(N)\ |A/(1)|; 2tt|/i 
1 / 1 



|A/(iV)| |A/(1) 



n) < 



n=l 



Allora 



1 - 

-Y 



2nihf{n) 



n=l 



< 



n\h\ \\Af(N)\ \Af( 



1 1 

+ 



1 \ N ^ 

m )+7r\h\^2mn)\ 



ir\h\ \N\Af(N)\ N\Af(l)\ 



N-l 



-™Eia/(»)|. 

71=1 

Tenendo conto delle ipotesi (|1.6|), dal Teorema 1.2.1 segue la tesi. □ 



Conseguenza immediata del teorema precedente (e, quindi, del criterio di 
Weyl) è il noto 

Teorema 1.2.4 (di Fejér). Sia f una funzione definita e differenziabile per 
x > 1 . Se f è monotona per ogni x > 1 ed inoltre 



lim f'(x) = 



x— >+oo 



:i-t) 



lim x\f'(x)\ = +oo 

X-++OC 

allora la successione (f(n)) ne n è u.d. mod 1. 
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Dimostrazione. Per il teorema del valor medio esiste r n e [n, n + 1] tale che 

A/(n) = /'(r») . 

Pertanto, (A/(n)) è monotona e soddisfa lim^oo A/(n) = in virtù dell'i- 
potesi (1.7). 



Quindi, per poter applicare il Teorema 1.2.3, rimane da verificare che A/(n) 
soddisfa la seconda in (1.6). Supponiamo dapprima che la /' sia decrescente. 
Allora, per n sufficientemente grande risulta che 

(n + l)f'(n + 1) - f\n + 1) = nf(n + 1) < nf(r n ) < nf{n). 



Passando al limite e tenendo conto delle ipotesi (1.7) e (1.8), otteniamo che 
lim n -+ 00 n\Af(rì)\ = +oo. 

Se, invece, /' è crescente, basta scambiare / con — /. 

In entrambi i casi possiamo applicare il teorema precedente ed ottenere la 
tesi. □ 

Esempi 1.2.5. Come abbiamo già osservato, il teorema di Fejér fornisce un 
criterio sufficiente per stabilire quando una successione di punti è u.d. mod 
1. Vediamone alcuni esempi proposti in [8]. 

1. (cm CT log T n) n > 2 , con a ^ 0, r 6 R, < a < 1. 

Vogliamo provare che la funzione f(x) = ax a log T x, con x > 2, associata 
alla successione di termine generale f(n) = an a log r n, verifica le ipotesi del 
teorema di Fejér. 

Osserviamo innanzitutto che si tratta di una funzione monotona. Infatti, si 
può provare che la derivata prima 



f'(x) = aax a log T x + 



ax u r log r 1 x 
x 



ax 



log r x(o~ logx + r) 



è sempre monotona. Inoltre, per x — > +oo questa derivata tende a zero e 
x\f(x)\ = \ax u log r_1 x(a\ogx + r)\ — > +oo . 



Pertanto, poiché la / verifica le condizioni (1.7) e (1.8), possiamo applicare 
il teorema di Fejér e dedurre che la successione è u.d. mod 1. 

2. (a log r n) n >i, con a^Oer > 1. 

Indichiamo con g(n) il termine generale della successione e poniamo g(x) = 
a log r x per x > 1. Allora 



olt log r 1 x 



g (x) 







(X — > +oo 



X 
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e inoltre si vede che g' è monotona e 

x\g'{x)\ = |a|rlog r_1 x — > +00 (x — > +00) . 



Per il Teorema 1.2.4, la successione (a log T n) n >i è u.d. mod 1. 
3. (cm log r n) n >2, con a 7^ e r < 0. 

Anche in questo caso possiamo verificare che sono soddisfatte le ipotesi 
del teorema di Fejér. Infatti, indicato con h(n) il termine generale della 
successione e posto h(x) = axlog T x per ogni x > 2, è facile osservare che 



, _ ars log a.-,. 
h (x) — a log x H = a log x{ log a; + r ) — )• (J (x — > +00) 

x 



che h! è monotona e che 

x\h'(x)\ = \a log T ~ 1 x(\ogx + r)\ — > +00 (x — >■ +00) . 
Quindi la successione data risulta u.d. mod 1. 



1.3 II caso multidimensionale 

Generalizziamo quanto detto finora al caso di successioni di vettori, partendo 
da alcune notazioni. 

Sia s un intero, s > 2. Siano a = (ai, . . . , a s ) e b = (61, . . . , b s ) due vettori di 
]R S . Diciamo che a < b (rispettivamente, a < b) se a, < òj (rispettivamente, 
< bi) per ogni % = 1, . . . , s. L'insieme dei punti x6R s tali che a < x < b 
(rispettivamente, a < x < b) sarà denotato da [a, b[ (rispettivamente, [a, b]) 
e verrà chiamato sottointervallo. 

Il cubo unitario s- dimensionale I s è l'intervallo [0, 1[, dove = (0, . . . , 0) e 
1 = (1,...,1). 

La parte intera di x = (xi, . . . , x s ) è = ( , • • • , [x s \ ) e la parte frazio- 
naria di x è {x} = ({xi}, . . . , {x s }). 
Sia (x n ), n — 1, . . . , una successione di vettori in W . 

Definizione 1.3.1. La successione (x n )„ g N è detta uniformemente distribui- 
ta modulo 1 {u.d. mod 1) m.W se 

j N s 
n=l j=l 

per tutti i sottointervalli [a, b[c I s . 
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Considerare il caso multidimensionale ci permette di fornire una defini- 
zione di uniforme distribuzione per successioni di numeri complessi. 

Definizione 1.3.2. Sia (z„) n eN una successione di numeri complessi. Siano 
Re z n = x n e Im z n = y n . Allora la successione (z n ) è detta uniformemente 
distribuita modulo 1 {u.d. mod 1) in C se la successione ((x n ,y n )) ne n è u.d. 
mod 1 in M 2 . 

Ricordando che il cubo unitario s-dimensionale chiuso I s è l'intervallo 
[0, 1], possiamo enunciare il seguente: 

Teorema 1.3.3. Una successione (x n ) n( z?q è u.d. mod 1 in M, s se e solo se 
per ogni funzione f continua a valori complessi definita su I s vale la seguente 
relazione: 



Infine, se per ogni x = (xi, . . . , x s ) e y = (yi, . . . , y s ) in 1R S denotiamo con 
(x, y) il prodotto scalare standard in W, ossia (x, y) = x%yi +£22/2 +■ • -+x s y s , 
possiamo scrivere il criterio di Weyl nel caso mutidimensionale, omettendone 
la dimostrazione che è essenzialmente simile a quella del caso s — 1. 

Teorema 1.3.4 (Criterio di Weyl). Una successione (x n ) n( z^ è u.d. mod 1 
in W se e solo se per ogni h G 7L S , h ^ 0, risulta 



1.4 Discrepanza 

Lo scopo di questo paragrafo è quello di presentare i concetti fin qui esposti 
da un punto di vista quantitativo, e non solo puramente qualitativo come 
fatto sinora. Per farlo introduciamo il concetto di discrepanza, termine pro- 
babilmente coniato da J. C. van der Corput, matematico olandese che per 
primo, nel 1935, pubblicò in [li] una procedura che consente di generare 
successioni u.d. mod 1 e che inoltre stabilì risultati quantitativi riguardanti 
il comportamento nella distribuzione di una successione. 
Questa quantità ci permette di distinguere tra "buone" e "cattive" succes- 
sioni u.d.. Infatti, tra le successioni u.d. mod 1 ve ne sono alcune che 
presentano una buona distribuzione, al contrario di altre che sono soltanto 
uniformemente distribuite. 




v->oo TV ^— ' 



lim ^^=0. 



n=l 
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Definizione 1.4.1. Sia x\, . . . , x^ una successione finita di numeri reali. La 
quantità 



D N = D N (x!, 



x N 



sup 

0<a<6<l 



1 N 



n=l 



;i.io) 



è detta discrepanza della successione. 

Per una successione infinita w di numeri reali (o per una successione di numeri 
reali finita contenente almeno N termini), la discrepanza Dn(w) si definisce 
come la discrepanza del segmento iniziale formato dai primi N termini di w. 



Teorema 1.4.2. La successione w è u.d. mod 1 se e solo se 

lim D N (w) = . 



Xil) 



Dimostrazione. Osserviamo che, per la Definizione 1.4.1 la (1.11) implica 
direttamente la (1.1), ossia la definizione di successione u.d. mod 1, per 
tutti gli intervalli [a, b[ C /. Quindi si tratta solo di dimostrare che ogni 
successione u.d. mod 1 soddisfa la (1.11). Per far ciò, scegliamo un intero 
positivo m, con m > 2. 

Per < k < m — 1 denotiamo con K l'intervallo K = I — , [■ 

Dalla Definizione 1.1.1 sappiamo che se una successione è u.d. mod 1 allora 

esiste iVo = N (m) tale che per ogni iV > Nq si ha 

- f 1 - -ì < 4 E < - f 1 + - 

v 7 n=l v 

Consideriamo, ora, un arbitrario intervallo J = [a, b[ C /. Chiaramente 
esistono due intervalli Ji, J2, unioni finite di intervalli J^, tali che J\ C J C J2 
e per cui sono verificate le seguenti relazioni: 

A(J)-A(J0< 2 



m 



A(J 2 ) - A(J) < - , 
m 

dove con A indichiamo la misura di Lebesgue unidimensionale. 
Allora, per ogni N > Nq risulta 



;i.i2) 
;i.i3) 



A(Ji 



/ 1 \ 1 1 

^ 7 n=l n=l 

1 J! / 1 \ 

< ^E^ 2 (K»<A(J 2 ) 1 + - 

n=l ^ 7 
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Da questa catena di disuguaglianze e dalle relazioni (1.12) e (1.13) si ottiene 
che 

A(J) - -) ( 1 - -) < ^ E Xj( W) < f A(J) + - Vi + - 

Infine, poiché J C /, dalla proprietà di monotonia della misura di Lebesgue 
segue che À(J) < 1 e quindi 



1 N 3 



n=l 



ViV > iV 



che in forma equivalente può essere scritta come 



1 N 



n=l 



3 2 
< - + — 



Poiché la stima precedente non dipende da J, possiamo considerare l'estremo 
superiore di ambo i membri su ogni possibile intervallo J C / ed ottenere 
così 



3 2 
D N (u) <- + 



m m 



2 ' 



Dal momento che la quantità — + -\ può essere resa arbitrariamente piccola, 
il teorema è completamente dimostrato. □ 

Osservazione 1.4.3. Vista l'arbitrarietà di J, il teorema precedente mette 
in luce un aspetto interessante sulle successioni u.d. mod 1, ossia che 



1 - 

lim J2xAM) = \(J) 



n=l 



uniformemente rispetto a tutti i sottointervalli J = [a, b[ di /. 

Dimostriamo, ora, un ulteriore teorema che mostra come Dn non possa 
convergere a zero troppo rapidamente. 

Teorema 1.4.4. Per ogni successione di finita X\,...,Xn di numeri reali 
risulta 



N 



< D N (xx,.. .,x N ) < 1 . 
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Dimostrazione. La disuguaglianza di destra segue in maniera evidente dalla 
definizione di discrepanza. Per dimostrare la prima disuguaglianza, fissiamo 
e > e scegliamo uno degli elementi della successione, che chiamiamo x per 
comodità. Consideriamo l'intervallo J — [x,x + e[f] I. Poiché x G J, si ha 



1 N 11 

^J2xAM)-x(j)>--x(j)>- 



e . 



n=l 



Dall'arbitrarietà di e segue la disuguaglianza di sinistra e ciò conclude la 
dimostrazione. □ 

A volte è utile restringere la famiglia di intervalli su cui calcoliamo l'e- 
stremo superiore nella definizione di discrepanza. Il tipo di restrizione più in- 
teressante consiste nel considerare intervalli della forma [0,c[, con < c < 1. 
La Definizione 1.4.1 viene modificata di conseguenza ed estesa nella maniera 
che segue. 



Definizione 1.4.5. Data una successione finita xi,x%, 
reali, definiamo star- discrepanza la quantità 



,%n di numeri 



D* 



D* N (xi, 



.xn 



sup 

0<c<l 



1 N 



n=l 



Analogamente a quanto succede per la discrepanza, la star-discrepanza 
D* N (u) di una successione cu infinita (o di una successione finita di almeno iV 
termini) si definisce come la star-discrepanza dei primi N termini di u. 

Il prossimo risultato ci permette di stabilire una connessione tra e D* N . 
Teorema 1.4.6. Tra le discrepanze Dn e D* N esiste la relazione 



D* N <D N < 2D* N 



;i.i4) 



Dimostrazione. La prima disuguaglianza segue direttamente dalle Definizioni 
ITXTl e [TA5l 

Per dimostrare la seconda, osserviamo che 



N 



N 



N 



X[a,b[({ X n}) = X[0,b[({ X nì) ~ Y X[0,a[(i X n}) (0 < a < 6 < 1) , 



n=l 



n=l 



n=l 



da cui segue che 



1 N 



71=1 



< 
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1 N 1 N 



n=\ 



n=l 



Passando all'estremo superiore otteniamo la ( |1.14 ). 
Corollario 1.4.7. La successione w è u.d. mod 1 se e solo se 



□ 



lim D* N {u 







Dimostrazione. E una conseguenza immediata di (1.11) e (1.14). 



□ 



Concludiamo questa sezione con una definizione che rende bene l'idea 
del carattere quantitativo della discrepanza. Appare chiaro dalle definizioni 
precedenti, infatti, che una successione di numeri reali è tanto più uniforme- 
mente distribuita quanto più velocemente converge a zero la sua discrepanza. 
E stato dimostrato che le successioni u di punti aventi la discrepanza asin- 
toticamente migliore sono quelle per cui D N (oo) è dell'ordine di ^j^- quando 
N — > oo. Vale, infatti, la seguente 

Definizione 1.4.8. Una successione u di numeri reali è detta successione di 
punti a bassa discrepanza se 



D N (u) = O (log N/N) {N ->■ oo) . 



Una successione speciale 

Esibiamo ora una successione molto particolare costruita sull'intervallo uni- 
tario: la successione di van der Corput. Ciò che la rende un interessante 
oggetto di studio è che essa è una successione infinita con bassa discrepanza. 
La costruzione di tale successione si basa sulla rappresentazione in base 2 dei 
numeri naturali e sulla funzione radice inversa, di cui ci accingiamo a fornire 
la definizione. 

Nel Capitolo 2 vedremo in dettaglio la rappresentazione di un numero natu- 
rale in una qualsiasi base b e la corrispondente funzione radice inversa. 

Definizione 1.4.9. Sia n EN. Se ri si scrive come YltLo a & \ dove a« 6 {0, 1} 
e M = \}og 2 n\ , la rappresentazione in base 2 di n è la sequenza di cifre 
[n] 2 = aMdM-i ■ ■ ■ Oo- La funzione 02 : N — > [0, 1], definita come 

M 
i=0 

è detta funzione radice inversa, e la sua rappresentazione in base 2 è 
[0 2 (n)] 2 = 0.a . . • a M -ia M - 
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Descriviamo brevemente questa funzione attraverso un facile esempio. 

Esempio 1.4.10. Consideriamo la rappresentazione binaria del numero na- 
turale 7. Si verifica facilmente che 

2 



i=0 

dove gli cij sono tutti uguali ale quindi [7] 2 = 111. Allora il valore della 
funzione radice inversa in 7 è 



la cui espressione binaria è [</>2(7)] 2 = 0.111. 

Ora, poiché il valore della funzione radice inversa può essere determinato 
per ogni numero naturale, risulta immediato considerare successioni numeri- 
che il cui termine generale sia proprio fain). 
La successione {(f>2 (n)} si ottiene nel modo seguente: 
1 . Si considera la rappresentazione binaria dei numeri naturali 



2. Si considera la funzione radice inversa di ciascuna rappresentazione 



Possiamo dunque presentare la seguente definizione. 

Definizione 1.4.11. La successione di van der Corput è definita come la 
successione di termine generale </> 2 (n), dove 02 è la funzione radice inversa. 

Enunciamo ora il risultato più importante presentato dallo stesso van der 
Corput su tale successione. 

Teorema 1.4.12. La discrepanza Dn(uj) della successione di van der Corput 
u = {02(^)}neN soddisfa 





i=0 



01,10,11,100,101,110,111,... 



113 15 3 7 
2'4'4'8'8'8'8'" 
o, equivalentemente, la loro espressione binaria 



0.1, 0.01, 0.11, 0.001, 0.101, 0.011, 0.111, 





per ogni N G N. 
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La conseguenza immediata di tale teorema è il seguente 
Corollario 1.4.13. La successione di van der Corput {4> 2 (n)} n( z?$ è u.d.. 



Dimostrazione. Segue immediatamente dal Teorema 1.4.2 perchè 1 °s(| v+1 ) è 



infinitesimo per iV — > oo. 



□ 



Ci sono varie generalizzazioni della successione di van der Corput. La 
prima estensione della procedura al caso bidimensionale avvenne per opera 
dello stesso van der Corput. Successivamente Hammersley [6] e Halton [5] 
hanno proposto altre generalizzazioni al caso di un ipercubo in qualsiasi 
dimensione, ma di questo ci occuperemo in maniera più estensiva nel terzo 
capitolo. 



Discrepanza multidimensionale 



Le definizioni di Dn e D* N possono essere estese in maniera piuttosto ovvia 
a successioni nello spazio s-dimensionale IR S . 



Definizione 1.4.14. Sia x 1; . . . , x^v una successione finita in 
panza e la star- discrepanza D* N sono definite come segue: 



D 



N 



LV(xi, . . . ,xjv) 



sup 

J 



1 N 

^5>j({x»})-a(j) 



n=l 



La discre- 



D 



N 



D 



sup 

J* 



1 N 

-^X^({x„})-A(J* 



n=l 



dove J varia tra tutti i possibili sottointervalli di I s della forma 
{x G M s : a < x < b} e J* varia tra tutti i sottointervalli della forma 
{x e l s : < x < b}, e A denota la misura di Lebesgue s-dimensionale. 

Infine, se la successione (x n ) è infinita, la discrepanza è quella del seg- 
mento iniziale formato dai primi N termini della successione, come abbiamo 
già visto nel caso unidimensionale. 

Vediamo come la relazione tra Dn e D* N in più dimensioni sia molto simile 
alla (1.14) che vale per s — 1. 

L'idea è quella di partire da uno dei sottointervalli J di I s e rappresentarlo in 
termini di intervalli del tipo J*. Descriviamo il procedimento per s = 2 perchè 
la generalizzazione a dimensioni superiori si ottiene in maniera analoga. 
Sia J il seguente sottointervallo di I 2 : 



J = {(xx,x 2 ) e I 2 : «i < xi < Pi A a 2 < x 2 < f3 2 } = [a x , 0i[x [a 2 , (3 2 [ , 
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con < a» < fa < 1 per i = 1,2. 
Allora 

J = {([0,/3 1 [x[0,/3 2 [)\([0,a 1 [x[0,/3 2 [)} 

\ {([0,^[x[0,a 2 [) \ ([0, Ol [x [0,05,0} = ( J i \ J 2) \ (J*3 \ Jl) ■ 

Di conseguenza, per le proprietà di cui gode la misura di Lebesgue, risulta 
A(J) = A(JD - A(J 2 *) - A(JJ) + A(J 4 *) 



A? 



N 



N 



N 



N 



n=l 



n=l 



71=1 



n=l 



n=l 



Con lo stesso ragionamento utilizzato nella dimostrazione del Teorema 1.4.6 
si conclude che 



D* N < D N < AD* N . 



;i.l6) 



Si può dimostrare che nel caso s-dimensionale la relazione tra e D* N 
diventa 



D* N <D N < 2 S D 



N 



N 



;i.i7) 



Analogamente a quanto fatto nel caso unidimensionale (Teorema 1.4.2 e Co 



rollano 1.4.7) si può dimostrare che una successione w è u.d. mod 1 se e solo 
se lim^oo D* N (w) = lim^oo D N (w) = 0. 



Infine, generalizzando il Teorema |1 .4.4 , ritroviamo anche nel caso multidi- 



mensionale la stima > jj- 
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Capitolo 2 

LiS-successioni di partizioni e di 
punti in [0, 1 



In questo capitolo prenderemo in esame dapprima una famiglia di successio- 
ni di partizioni nell'intervallo unitario, introdotte per la prima volta da S. 
Kakutani in [7], ed il loro rapporto con la teoria dell'uniforme distribuzio- 
ne. Successivamente, focalizzeremo la nostra attenzione sulle LS'-successioni, 
introdotte per la prima volta da I. Carbone in [2J. Seguendo la linea di 
questo lavoro, prenderemo in esame le successioni di punti associate alle LS- 
successioni di partizioni, vedremo quale relazione intercorre con la successione 
di van der Corput, di cui abbiamo già studiato la costruzione nel Paragrafo 



1.4| e di cui sappiamo che la discrepanza 'è ottimale ed, infine, introdurre- 
mo un algoritmo "veloce" [3] completando il capitolo con la presentazione 
di un programma da noi implementato. L'utilizzo dell'aggettivo veloce non 
è casuale, ma nasce dall'esigenza di confronto tra quest'ultimo e l'algoritmo 
presentato in [2J, rispetto al quale ha una complessità computazionale meno 
elevata. 



2.1 Successioni di partizioni e raffinamenti 

Definizione 2.1.1. Una partizione n dell'intervallo I = [0, 1[ è una famiglia 
di suoi sottointervalli a due a due disgiunti, determinati da una successione 
crescente di punti {yo, yi, ■ ■ ■ ,yk}, ossia 

= {[yi,y l+ i[: 0<i<k-l} 

con I = U k i Zo[y i ,y l+ i[ ■ 
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Il concetto di uniforme distribuzione si estende in maniera naturale alle 
successioni di partizioni dell'intervallo unitario nel modo seguente. 

Definizione 2.1.2. Una successione di partizioni {7r n } neN di [0, 1[, dove n n = 
{[vt^iVÌ+iV- 1 < z < k(n)}, è detta uniformemente distribuita (u.d.) se 

per ogni funzione continua / definita su [0, 1]. 

Come abbiamo già visto nel Capitolo 1, la definizione precedente può 
anche essere espressa nella seguente forma equivalente. 

Definizione 2.1.3. Diciamo che la successione di partizioni {vr Tt } n6 N con 
= { [Ui^ , Vi+i [ '■ 1 < i < k(n)} è u.d. se per ogni coppia di numeri reali 
a, b, con < a < b < 1, si ha 



k(n) 

k(n) ^ 



lim jLT,X[aM(y < i n) ) = b-a (2.2) 

n—>-co ■ 



o, equivalentemente, 

lim t^-t X[o,c[(^ (n) ) = c per ogni < c < 1 . (2.3) 



k(n) 



n^oo k(n) 

v ' 1=1 



Ora descriviamo brevemente la procedura introdotta da Kakutani nel 
1976 per generare una classe particolare di successioni di partizioni u.d. 
dell'intervallo [0,1]. 

Definizione 2.1.4 (a-raffinamenti di Kakutani). Sia a G ]0, 1[ e sia ti = 

{[ìJiiUì+A '■ 1 < i < A; — 1} una partizione di [0, 1[. Chiamiamo a -raffinamento 
di ir, e lo indicheremo con air, la suddivisione in due parti dei soli intervalli 
di lunghezza massima di tt, proporzionalmente ad a e 1 — a, rispettivamente. 
Se denotiamo con a n 7r l'a-raffinamento di a n ~ l ii per ogni n G N, la succes- 
sione {/tri} n6 N ottenuta per successivi a-raffinamenti della partizione banale 
u> — {[0, 1[} è detta successione di partizioni di Kakutani. 

Osserviamo che per a — |, la successione di partizioni di Kakutani coin- 
cide con la partizione binaria. 
Fu lo stesso Kakutani a dimostrare che 
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Teorema 2.1.5. La successione di partizioni di Kakutani {n n } n( z N è u.d. per 
ogni a G ]0, 1[. 

Al fine di illustrare le LS'-successioni di partizioni, e le successioni di 
punti ad esse associate, è opportuno presentare il concetto di p-raffinamento 
introdotto da A. Volcic in [12], che generalizza quello di a-raffinamento di 
Kakutani. 

Definizione 2.1.6 (p-raffinamento). Sia p una partizione non banale di [0, 1[. 
Definiamo p-raffinamento di una partizione ir di [0, 1[, e lo indicheremo con 
p7r, la suddivisione di tutti gli intervalli di ir che hanno lunghezza massima 
in maniera omotetica rispetto a p. 

Osservazione 2.1.7. Se consideriamo la partizione p = {[0, a[, [a, 1[} e pren- 
diamo 7T = u), allora il suo p-raffinamento coincide con l'a-raffinamento di 
Kakutani. 

Inoltre, come nel caso della procedura di Kakutani, possiamo considerare il 
p-raffinamento di pn ed indicarlo con p 2 7r e così via, indicando con p n ir il 
p-raffinamento di p n ~ Xr n . 

Definizione 2.1.8. Data una partizione non banale p di [0, 1[, la successione 
dei p-raffinamenti {p""n:}nm della partizione 7r si definisce come la successione 
di partizioni ottenuta mediante i successivi p-raffinamenti di ir. 

Se con {p™a;} ng N indichiamo la successione di successivi p-raffinamenti 
della partizione banale u, vale il seguente importante risultato (|12j. Theorem 
2.7). 

Teorema 2.1.9. La successione {p n u} ne fq è u.d.. 

Vediamo infine come è definita la discrepanza per una successione di 
partizioni e quando questa può dirsi ottimale. 

Definizione 2.1.10. La discrepanza di una successione di partizioni {7T n } ne pj, 
con 7r n = { [s/i j 2/ì+i [ : 1 < i ^ k(n)}, è definita come 

k(n) 

^X[a,6[(vj n) ) - ( b ~ a ) ■ 

i=i 

È stato dimostrato che le successioni di partizioni {7r n } ng N aventi la di- 
screpanza asintoticamente migliore sono quelle per cui D(7r n ) "è dell'ordine 
di quando n — y oo. 

È stata introdotta, perciò, la seguente 



D(7T n ) = sup 

0<a<ò<l 
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Definizione 2.1.11. Una successione di partizioni {vr„} ngN , dove n n = 
{[Ui^ìUi+A'- 1 < i < k(n)}, è detta a òassa discrepanza se 

D(n n ) — 0(1/ k(n)) (n-»oo). 

2.2 L/S-successioni di partizioni 

In questo paragrafo prenderemo in esame una classe particolare di p-raffina- 
menti della partizione banale u dell'intervallo [0, 1[ : le LS'-successioni di 
partizioni, che costituiscono una classe numerabile di successioni, al variare 
di L ed S in N. 

Definizione 2.2.1. Fissiamo due interi positivi L ed S e sia /3 il numero reale 
positivo che verifica l'equazione + S(5 2 = 1. Indichiamo con p^g la parti- 
zione definita da L intervalli lunghi aventi lunghezza f3 e da S intervalli corti 
aventi lunghezza (3 2 . Chiameremo LS -successione di partizioni la successione 
{pls^ìneN (brevemente, {p^sìnen) dei successivi p Lj s-raffinamenti della par- 
tizione banale u> dove, ovviamente, con p\ s denotiamo il p^s-raffinamento 

di pXI 

Illustriamo meglio questa costruzione attraverso alcuni esempi. 

Esempi 2.2.2. 1. Consideriamo L — S — 1 e la corrispondente 1, 1- succes- 
sione di partizioni {Pii\ n e^ determinata dalla relazione j3 + (3 2 = 1, la cui 
soluzione in [0, 1[ è il numero ci/urco j3 — - . Lei prima partizione consterà 
semplicemente di due intervalli di lunghezza (3 e (3 2 , rispettivamente, ossia 

P \ tl = {[o,fi[, \p,i\}, 

come mostra il seguente grafico: 



Per costruire la seconda partizione, si suddivide l'intervallo di lunghezza 
maggiore (il primo) proporzionalmente a (3 e /3 2 in modo da ottenere da esso 
un intervallo lungo ed uno corto, e cioè 



p lil = {[0,/3 2 [, [/3 2 ,/3[, [/?,![} 



2.2 LS-successioni di partizioni 



24 



/3 2 P 1 



A questo punto, come si può facilmente osservare, ci sono due intervalli 
di lunghezza massima, il primo ed il terzo, per cui si procede alla loro sola 
suddivisione. Otteniamo così la terza partizione 

PÌ,i = {[0, P\ [P\ P% [P\ P[, [P, P + P% [P + P\ 1[} , 
la cui rappresentazione grafica è 



o p 3 p 2 P p + p 3 1 



e così via in modo da ottenere le altre. 
Sottolineiamo che questo è un esempio di /3-raffinamento alla Kakutani, cioè 
{Pi,i}nen è la successione di partizioni di Kakutani ottenuta per successivi 
/3-rafhnamenti di u. 

Passiamo ora ad un esempio un po' più complicato. 

2. Consideriamo la successione di partizioni {p^iWn corrispondente a 
2/3 + P 2 = 1. Al primo passo avremo una suddivisione dell'intervallo uni- 
tario in due intervalli lunghi ed uno corto, in modo da ottenere la seguente 
partizione 

pì,i = W,P[, [PM [2p,i[} , 



p 2(3 1 

Per ottenere la seconda partizione si dovranno pertanto dividere i primi 
due intervalli, in modo da ottenere da ciascuno di essi due intervalli lunghi e 
uno corto. Si avrà dunque 

PÌ,i = {[0,P 2 1[P 2 ,2P 2 1[2P 2 ,PI[P,P + P 2 [, 

\fi + P*,0 + 2P 2 [, [P + 2P 2 ,2P[, [2/5, 1[} , 



che è graficamente rappresentato qui di seguito. 
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/3 2 2/3 2 P (3 + f3 2 (3 + 2^2/3 1 



Al passo successivo gli intervalli da dividere saranno cinque, ciascuno di 
lunghezza P 2 . 

Iterando questo processo otteniamo le partizioni successive. 

3. Infine, presentiamo la successione di partizioni {p" 2 }neN relativa a 
P + 2(3 2 = 1. Sebbene (come vedremo in seguito) a differenza delle due pre- 
cedenti non si tratti di una successione a bassa discrepanza, tale esempio sarà 
utile per la comprensione dell'algoritmo nella sua formulazione generale. 
La prima partizione si ottiene dividendo l'intervallo unitario in un intervallo 
lungo di lunghezza fi — \ e in due intervalli corti di lunghezza (3 2 = |, nel 
modo seguente: 



o P p + p 2 1 



cioì 

PÌ,2 = W,P[, + A [/3 + /3 2 ,l[}. 

Al secondo passo dobbiamo suddividere l'unico intervallo di lunghezza mas- 
sima (3 proporzionalmente a (3 e /3 2 , ottenendo così 

ph = ito, n [p\ p 2 + p% [p 2 + p\ pi ip, p+p%[p+ p 2 , i[> , 

che è rappresentato così: 



^2 (3 2 + (3 3 p /3 + /3 2 1 



Procedendo in maniera analoga a quanto fatto finora si ottengono tutte 
le altre partizioni. 
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Prendiamo in considerazione alcune proprietà di tali successioni, che met- 
tano in luce il loro carattere ricorsivo e vediamo come tra esse meriti una par- 
ticolare attenzione la successione ottenuta considerando la coppia di interi 
L = S = 1. 

Osservazioni 2.2.3. (i) Notiamo innanzitutto che per ogni partizione p\ s 
contiene solo due tipi di intervalli: quelli lunghi e quelli corti, i primi di misura 
/3 n , i secondi /3 n+1 . Se con t n indichiamo il numero totale di intervalli di p\ s , 
allora vale la seguente equazione alle differenze finite, lineare, omogenea, del 
secondo ordine 

t n = Lt n -x + St n - 2 , (2.4) 
con condizioni iniziali to — 1 e ti — L + S. 

Se indichiamo con l n il numero di intervalli lunghi di s e con s n il numero 
dei suoi intervalli corti, valgono le seguenti relazioni: 

t n In ~\~ S n , l n Ll n —\ ~\- S n —i , S n Sl n —\ . 

Infatti l'n-esima partizione si ottiene dalla somma dei suoi intervalli lunghi e 
dei suoi intervalli corti. In particolare, gli intervalli lunghi della partizione ri- 
esima sono dati dalla suddivisione degli l n -\ intervalli lunghi della partizione 
Pls ' proporzionalmente in L intervalli lunghi, sommati agli intervalli corti 
della partizione precedente p^s cne diventano lunghi al passo successivo, 
mentre gli intervalli corti sono quegli intervalli ottenuti suddividendo gli Z n _i 
intervalli lunghi della partizione p^g, proporzionalmente in S intervalli corti. 
Allora, utilizzando queste relazioni si ha che 

t n l n -\- s n Ll n —i -\- s n —i -\- Sl n —\ 

= Un-\ + Sl n -2 + SXZ n „2 + Ss n -2 — 
— L{l n _i + + S(l ri ^2 + ■S„_2) = 

(ii) Dalla (2.4) segue immediatamente che la successione {tn^rigN) riel caso 
L — S — 1, è la ben nota successione di Fibonacci. D'altra parte abbiamo già 



evidenziato nell'Esempio 2.2.2 che la successione {pii} n eN è una successione 
di Kakutani. Per questi due motivi insieme, la 1, 1-successione di partizioni 
corrispondente è stata chiamata in [2] successione di partizioni di Kakutani- 
Fibonacci. 

In ([2], Theorem 2.2) è stato dimostrato, altresì, che per L > S tutte le 
LS'-successioni di partizioni hanno discrepanza ottimale, e cioè D(p2 s ) = 

O ( j- ) per n — > oo, ed è stata data una stima per la discrepanza in tutti i 
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casi. Osservato che L > S equivale a S < L+l, e ricordato che L(3 + S(3 = 1 
si ha il seguente 



Teorema 2.2.4. (i) Se S < L + 1 esistono due costanti positive C\ e ci tali 
che per ogni nGN 

f < D(pls) < f . 
(ii) Se S = L + l esistono due costanti positive c 3 e c 4 tali che per ogni n e PJ 

log*™ n , n . logt n 

C 3 — < £>(Pl,s) < C 4-7 • 

(ii) Se S > L + l esistono due costanti positive c 5 e c 6 ta/z c/ie per ogni n e PJ 
dove 7 = 1 + < l 

' log/3 

Sottolineiamo, infine, che grazie a questo teorema, per la prima volta in 
letteratura è stata presentata una stima per la discrepanza di una successione 
di Kakutani (che risulta essere, inoltre, a bassa discrepanza). 



2.3 L/S-successioni di punti 

Passiamo ora allo studio delle successioni di punti (anche'esse introdotte in 
[2]) associate alle L^-successioni di partizioni. Lo scopo è quello di presen- 
tare il riordinamento "alla van der Corput" della successione di punti asso- 
ciata a ciascuna LS'-successione di partizioni, in modo che la corrispondente 
successione di punti abbia discrepanza migliore possibile. Con la locuzione 
"alla van der Corput" intendiamo il processo di riordinamento che ricorda la 



costruzione della successione di van der Corput, studiata nel Paragrafo 1.4 

Definizione 2.3.1. La LS -successione di punti {££s}neN associata alla suc- 
cessione di partizioni {p£ s}neN è definita come segue. I primi ti punti so- 
no esattamente gli estremi sinistri degli intervalli di p\ s , presi nell'ordi- 
ne lessicografico. Tale insieme ordinato sarà denotato con s ed i suoi 

punti saranno indicati con Ci )---)Cti ■ m generale, per n > 1, se con 

A2 S = Kj n \ • • • ' indichiamo i t n punti della partizione p'1 s , scritti 

nell'ordine giusto, i punti della partizione successiva p 7 ^ si ottengono appli- 
cando agli estremi sinistri degli intervalli lunghi della partizione precedente 
le seguenti famiglie di funzioni 

^ +1 \x) =x + i(3 n+1 e ^ 1] (x) =x + L(3 n+1 + j(3 n+2 , (2.5) 
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con 1 < z < L e 1 < j < S - 1. 

Precisamente, i t n+ \ punti della partizione p 1 ^ sono ordinati nel modo 
seguente: 



n+l _ / An) An) 
L,S ~ \yi > ■ ■ ■ J 'Sin ' 

(n+l)/ c .(n)N (n+l)/ c (n)N (>i+l)/ t (n)\ (n+l)/ t (n^ 



, n ( n + 1 )(A n )\ ,A n+L > < C^ n >\ ,A n+1 > I A n >\ ,rS n+1 > ( A n >\ 

Va Isi J>--->Va \^l n )ì---:YL Isi )ì---iYL \H n ) 

ir ,( n+1 )(A n )\ (n+l),An)s , ( n+1 ) (A n )\ ,o( n+1 ) (c( n h\ 

YL,l 1^1 )ì---ìYL,1 lsz n J) • • • j <Pl,S-1\$1 li ■ ■ ■ i YL,S-l\Sl n )) ■ 

In ([2J, Theorem 3.9) è stato dimostrato che D(^ s , . . . , £g fS ) = O (^jf) 
per n — » oo, se L > S, e cioè, che a successioni di partizioni {p£ sì a bassa 
discrepanza corrispondono successioni di punti {££ s } a bassa discrepanza. 
Richiamiamo, a questo punto, il seguente 

Teorema 2.3.2. (i) Se S < L + 1 esiste una costante positiva ki tale che 
per ogni ÌVgN 

D N (e L ,s,e L , s ,---^Zs)<h 1 ^. 

(ii) Se S = L + 1 esiste una costante positiva ki tale che per ogni iVeN 
n iti p2 pN s < u lo ^ N 

(ii) Se S > L + 1 esiste una costante positiva k 3 tale che per ogni N e N 

logiV 



dove t = 1 + ^ < 1 



log/ 



Per maggiore chiarezza, facciamo vedere come si ottengono le successioni 
di punti corrispondenti alle successioni di partizioni presentate negli Esempi 

[2X2 

Esempi 2.3.3. 1. Analizziamo dapprima la successione di punti {£i i}neN- 
L'insieme di punti corrispondente alla partizione è costituito esattamente 
dagli estremi sinistri dei due intervalli [0,(3[, [(3, 1[, cioè 



Ai,i= (0,0) = • 



A questo punto, osserviamo che in questo esempio non ci sono funzioni (2.5) 
del tipo <p£~t- in quanto abbiamo un solo intervallo corto e quindi j = 0, ma 
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solo la funzione <p™ . Poiché li — 1, applicando la funzione (p\ al punto 
zero si ha 

Al = (0,/3,^ 2) (0)) = (0,/3,/3 2 ) = ■ 

\ (3) 

Poiché Z 2 = 2, applichiamo la funzione (p\ ai punti e /3, e precisamente 

Ali = (0,/3,/3 2 ,^ 3) (0),^ 3) (/3)) 

= (0,/3,/3 2 ,/3 3 ,/3 + /3 3 ) = • 

Applicando poi la funzione (fi ai primi tre punti (perchè I3 — 3), otteniamo 
Af j di cui riportiamo il grafico che segue. 

1° 6° 4° 3° 8° 2° 7° 5° 

I 1 1 1 1 1 1 1 1 

/3 4 /3 3 P 2 /3 + /3 3 P P + /3 4 (3 + p 3 

Figura 2.1: Primi 8 punti di {C^ilneN 

Così continuando otteniamo la successione di punti {£"i} n eN associata 
alla successione di punti di partizioni {p" x} n eN- 

Questa successione di punti è stata chiamata in [2] successione di punti 
Kakutani- Fibonacci, dal momento che discende dalla successione di partizioni 
di Kakutani-Fibonacci. 

2. Consideriamo, adesso, la successione di punti {£2 i}n.eN- Anche in questo 
caso non ci sono funzioni del tipo P^l^ da considerare, ma soltanto le fun- 
zioni del tipo (p\ n+1 \ con i — 1, 2, da applicare ai primi l n punti dell'insieme 
A% v II primo insieme di punti è 

= (0, 0, 2/3) = . 

Poi, essendo li = 2, se applichiamo le due funzioni e (p^ a i primi due 
punti di K\ 1 abbiamo 

A|i = (0,/3,^S 2) (0),^ 2) (/3),^ 2) (0),^ 2) (/3)) 
= (o,/3,2/3,/3 2 ,/3 + /3 2 ,2/3 2 ,/3 + 2/3 2 ) 

_ (fi ti t3 t4 ^5 t6 {7 \ 
— \ S2,l5 S2,l> S2,l! S2,l) S2,l) S2,l' S2,l ) ■ 



2.3 LS-successioni di punti 



30 



Calcoliamo ancora un insieme di punti, quello relativo alla partizione p 21 . 

Applicando le funzioni tp^ e tp^ ai primi cinque punti di A 2jl (visto che 
I2 = 5), si ha 

= (o,/3,2/3,/3 2 ,/3 + /3 2 ,2/3 2 ,/3 + 2/3 2 , 

pf (0) , <pf {fi) , <pV> (2/3) , <pf (/3 2 ) , <pf {fi + /3 2 )) 

= (o,/3,2/3,/3 2 ,/3 + /3 2 ,2/3 2 ,/3 + 2/3 2 , 

/3 3 , f3 + /3 3 , 2/3 + /3 3 , /3 2 + /3 3 , /3 + /3 2 + /3 3 , 

2/3 3 , /3 + 2/3 3 , 2/3 + 2/3 3 , /3 2 + 2/3 3 , (3 + /3 2 + 2/3 3 ) 

_ /Vi £-2 e 3 e 4 t 5 c 6 e 7 e 8 
~~ l S2,l) S2,l) S2,l) S2,l) S2,l) S2,l) S2,l) S2,l) 

_ e 9 tf 10 eli £12 c 13 tf 14 t 15 t 16 t 17\ 
— S2,l) S2,l; S2,l; S2,l) S2,l) S2,l> S2,l) S2,l) S2,l J • 

Tali punti hanno la seguente rappresentazione grafica: 



8° 13°4° 



ll c 



16°6° 



14°5° 



12° 17°7° 



H ; 

3 + /3 2 
2/3 3 /3 + /3' 



+ 



10° 15° 



H — I 



,3 2/3 3 /3 2 /3 2 + /3 3 I 2/3"= 
■ 2/3 3 



3 2 + p 3 2/3 , 
/3 + 2/3 2 2/3 - 



Figura 2.2: Primi 17 punti di {£ 2) i} n eN 

Iterando questo processo riusciamo a determinare tutti i punti della suc- 
cessione {^2,i}neN- 

3. Consideriamo infine l'esempio delle successioni di punti {£™ 2 }neN e de- 
terminiamo l'insieme di punti associato a ciascuna delle tre partizioni p\ 2 , p\ 2 
e p\ 2 - Il primo insieme di punti si determina facilmente sempre considerando 
gli estremi sinistri degli intervalli di p\ 2 , cioè 

A{ 2 = (0, /3, /3 + /3 2 ) = (£ì, 2 ,£i 2 2,£i 3 2 ) • 

Poiché /1 = 1, dobbiamo applicare le funzioni p^ e pf^ solo al primo punto, 
e cioè allo zero. Per cui 

A? >2 = (0,/3,/3 + /3 2 ,^ 2) (0),^ 2 ;(0)) 
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= (o,/3,/3 + /3 2 ,/3 2 ,/3 2 + /3 3 ) 

= ^£l,2! £l,2>£l,2> £l,2> £1,2) • 
Poiché l 2 = 3, procedendo come nei casi precedenti si trova che 

A?, 2 = (o,/3,/3 + /3 2 ,/3 2 ,/3 2 + /3 3 , 

^(0), ^S 3) (/3),^S 3) (/3 + /3 2 ), </>ì>), ^S 3 ì(/3 + /? 2 )) 

= (o,/3,/3 + /3 2 ,/3 2 ,/3 2 + /3 3 ,/3 3 ,/3 + /3 3 , 

/3 + /3 2 + /3 3 , /3 3 + f3\ /5 + /3 3 + /3 4 , /3 + /3 2 + /3 3 + /3 4 ) 

_ /ti £-2 t3 t4 t5 t6 t7 t8 £-9 tlO tll A 
— l Sl,2) Sl,2> Sl,2) ?1,2) ?1,2) £l,2> Sl,2> Sl,2) Sl,2) M.,2) Sl,2 ) ■ 

E cos via. 

L'insieme di punti A 3 2 si pu rappresentare nel modo seguente: 



1° 


6° 9° 4° 


5° 


2° 


7° 10° 3° 8° lì 


L° 


1 




1 1 1 

03 pi + p4 02 


1 +/3 3 


P 


1 i 1 1 ] 

/3 + fS 3 \ I3 + /3 2 /3 + /3 2 +/3 3 


; 1 



+ /3 3 + /3 4 + /3 2 + /3 3 + /3 4 



Figura 2.3: Primi 11 punti di {£" 2 }neN 

Nel prossimo paragrafo presenteremo un metodo meno complesso per il 
calcolo dei punti associati alle LS'-successioni di partizioni. 

2.4 Procedimento "alla van der Corput" 

In questo paragrafo costruiamo un certo numero di funzioni definite sull'in- 
tervallo unitario, dipendenti da L ed S, che ci permetteranno di determinare 
gli insiemi di punti A£ s e > di conseguenza, i punti di {££s}neN [3]. Per de- 
scrivere la loro costruzione abbiamo bisogno di introdurre la definizione di 
rappresentazione di un numero rispetto ad una base e una nuova famiglia di 
funzioni sull'intervallo unitario. 
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Rappresentazione dei numeri naturali 

Definizione 2.4.1. La rappresentazione di un numero è detta posizionale 
se esso è rappresentato da una sequenza di cifre il cui valore dipende dalla 
posizione all'interno della sequenza. 

Il nostro sistema di numerazione è posizionale decimale, ossia utilizza una 
rappresentazione in base 10, in cui il valore delle cifre, da a 9, è determinato 
dalla posizione che assume nella sequenza, per cui 1234 è diverso da 4321. 



A questo punto, data l'importanza che la rappresentazione in base b rico- 
prirà nel seguito della tesi, chiariamo come ottenere tale rappresentazione. 

Definizione 2.4.2 (Numeri in base b > 2). Siano b > 2 un intero ed n un 
numero naturale. Per rappresentare n in base b, si divide n per b e si può 
presentare uno dei tre casi: 

1. n = 

2. n < b 

3. n > b . 



Se n = 0, oppure n < b, abbiamo finito. Se, invece, n > b, siano ao e qo, 
rispettivamente, il resto ed il quoziente della divisione di n per b. Proseguia- 
mo dividendo qo per b e siano ai e q±, il resto ed il quoziente della divisione 
di qo per b, rispettivamente. Iteriamo questo procedimento fino a che non 
troviamo un resto a M nullo e, di conseguenza, una successione di resti del 
tipo aoai . . . om-iCìm- Allora la rappresentazione di n in base b è 

[n]b = aM^n-l ■ ■ ■ OiOq . 

Ovviamente, il numero n in base dieci è 



Al 



n = \n\ 1Q 



a M b M + aM-i^ 1 ' 1 + ... + ai b + a = J2 a i b * ■ ( 2 - 6 ) 



i=Q 



Tenendo conto della definizione precedente, possiamo ora estendere la 



Definizione 1.4.9 di funzione radice inversa rispetto ad una generica base b. 



Definizione 2.4.3. Siano b > 2 un intero ed n G N. Se la rappresentazione 
in base b di n è [n] b = om^m-i ■ ■ ■ a o, allora la funzione (fib '■ N — > [0,1], 
definita come 

M 
i=0 
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è detta funzione radice inversa, e la sua rappresentazione in base b è 
[4>b{n)]b = 0.a . . . ùM-idM- 

Esempio 2.4.4. L'espressione del numero 70 in base 3 è 2121. Infatti, 
70 : 3 = 23 con resto do = 1, 23 : 3 = 7 con resto 02 = 2, 7 : 3 = 2 
con resto a2 = le2:3 = con resto a 3 = 2. Il resto successivo è nullo e 
quindi il processo si arresta con 03. 
Pertanto [70] 3 = 2121 e 



[70] 



10 



3 

i=0 



1 • 3° + 2 ■ 3 1 + 1 • 3 2 + 2 • 3~ 



Inoltre 



. v — v CLj 12 1 



i=0 



2 

3 1 



e [0 3 (7O)] 3 = 0.1212. 



Sappiamo che la successione di van der Corput si ottiene considerando 
la rappresentazione in base 2 dei numeri naturali e invertendo l'ordine delle 



sue cifre (si veda Definizione 1.4.9 e Definizione 1.4.11). Se poniamo b = 2 e 



a = |, si 


hanno le 


seguenti relazioni: 










0^ 


[0] 2 = 





-> 0.0 = 


[0 2 (O)] 2 -► 







= 02 (0) 


1 


[1]2 


01 


->■ 0.10 = 


[0 2 (1)] 2 -► 


Q 




= 02(1) 


2 ->■ 


[2] 2 


10 


->■ 0.01 = 


[0 2 (2)] 2 -► 


a 2 




= 02 (2) 


3 -> 


[3]2 = 


11 


-> 0.11 = 


[<fe(3)] 2 -> 


a - 


Ha 2 


= 02 (3) 


4^ 


W2 = 


100 


-> 0.001 = 


[0 2 (4)] 2 -> 


a 3 




= 02(4) 


5 -> 


[5]2 = 


101 


-> 0.101 = 


[0 2 (5)] 2 -> 


a - 


h a 3 


= 02(5) 


6 -> 


[6]2 = 


110 


-> 0.011 = 


[0 2 (6)] 2 -> 


a 2 


+ a 3 


= 02(6) 


7^ 


[7] a = 


111 


-> 0.111 = 


[0 2 (7)] 2 -> 


a - 


h« 2 +a 3 


= 02(7) 


8 


[8]2 = 


1000 


->■ 0.0001 = 


[0 2 (8)] 2 -> 


a 4 




= 02(8) 


9 


[9] 2 


1001 


->■ 0.1001 = 


[0 2 (9)] 2 -> 


a - 


Ha 4 


= 02(9) 


IO-» 


[10]a = 


1010 


->■ 0.0101 = 


[0 2 (1O)] 2 -> 


a 2 


+ a 4 


= 02(10) 



È interessante e utile notare a questo punto che tale successione si ottiene 
riordinando i punti determinati dalla successione di partizioni di Kakutani 
per a = \ (o, altrimenti detto, la successione di partizioni binaria). 

È immediato osservare che le potenze di a dipendono dalla posizione decimale 
e questo ci permette di conoscere l'n-simo termine della successione di van 
der Corput utilizzando [n] 2 . 
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Questa continua ad essere un'osservazione valida anche se consideriamo basi 
diverse da 2. 

A questo punto ci chiediamo se esista una corrispondenza simile per i punti 
delle successioni {^LsIneN, e cioè tra le potenze di (3, dove L/3 + S(5 2 — 1, e 
la rappresentazione in qualche base. Scopriremo che in un certo senso è così. 
Partiamo dagli esempi già considerati nel paragrafo precedente, e vediamo 
come riottenere i punti delle successioni attraverso l'introduzione di nuove 
funzioni. 

Esempio 2.4.5. In questo esempio prendiamo in considerazione la succes- 
sione di punti di Kakutani-Fibonacci {£i\} n eN- Partiamo dalla successione 
di partizioni di Kakutani-Fibonacci {p^ijneN e definiamo le seguenti funzioni 



e 



V> :[0,1[— ► [0,/3[ 
x \-> j3x 



^:[0,/9[— ► 
x h-> f3x + (3 . 



(2.7) 
(2.8) 



Vedremo immediatamente che tutti i punti della successione {^"xjneN si ot- 
tengono attraverso composizioni successive di e ^. 

Osserviamo che le composizioni possibili tra le due funzioni sono ipo o i/j q , 
ipo o ifj 1 e ipi o i/jq, mentre ip\ o ip l non è consentita. Indicheremo con ip^ le 
composizioni ipi o ipj. Osserviamo inoltre che l'insieme dei punti relativi alla 
prima partizione è dato da 

Aì,i = (eì > i,e?,i) = (^(o)^i(o)) = (o,/3). 

Applicando queste funzioni (ovviamente quelle possibili) ai punti trovati e 
osservando che V'ìo(O) — V'i(O) P er °g n i * = 0? 1> s i na proprio 

Ali = = (Vo(0),Vi(0), Voi (0)) 

= (Voo(0),Vio(0),Voi(0)) 

= (0,/3,/3 2 ) • 

Iterando nuovamente questo procedimento otteniamo 

A?,! = (^,1, ^,1^,1,^,1^1,1) = (Vo(0),Vi(0),Voi(0),Vooi(0),Vioi(0)) 
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= (^ooo(O), ^ioo(O), Voio(O), Vooi(O), ^ioi(O) 
= (o,/3,/3 2 ,/3 3 ,/3 + /3 3 ) . 

Si vede chiaramente come la composizione ricorsiva delle funzioni definite in 
precedenza, che dipendono esclusivamente da L ed S, dia luogo alle successio- 
ni di punti associate alle LS- successioni di partizioni, ma la caratteristica sor- 
prendente consiste soprattutto nell'ordinamento dei punti, i quali presentano 
già un ordinamento "alla van der Corput". Questo appare particolarmente 
evidente se consideriamo l'insieme di punti successivo 

— ^14)^1,1^14)^1,1)^1,1)^1,1)^1,1)^1,1^ 
= (Vo (0) , Vi (0) , Voi (0) , Vooi (0) , Vioi (0) , 

Voooi (0) , Viooi (0) , Voioi (0)) 
= (Voooo(0), Viooo(O), Voioo(O), Vooio(O), 

Vioio (0) , Voooi (0) , Viooi (0) , Voioi (0)) 
= (0,/3,/3 2 ,/3 3 ,/3 + /3 3 ,/3 4 ,/3 + /3 4 ,/3 2 + /3 4 ) . 

Così continuando si ottiene la successione di punti {C^ijneN- 

A questo punto si tratta di osservare che esiste la seguente relazione tra le 

potenze di (3 e la rappresentazione binaria dei numeri naturali: 

[o] 2 - oo^ o.oo = [<M0)] 2 — > o 

1 — ► [1] 2 = 01 — > 0.10 = [<h0)h — > /3 

2 — > [2] 2 = 10 — > 0.01 = </> 2 (2)] 2 — -> /3 2 

3 — ► [3] 2 = 11 — >• 0.11 = [</> 2 (3)] 2 — ► nessun /3 

4 — > [4] 2 = 100 — > 0.01 = [<fe(4)] 2 — -> /3 3 
5^ [5] 2 = 101— > 0.101 = [<fe(5)] 2 — > /3 + /3 3 

6 — >• [6] 2 = HO — > 0.011 = [<^ 2 (6)] 2 — > nessun /3 

7 — > [7] 2 = IH — >• 0.111 = [^ 2 (7)] 2 — ► nessun /3 

8 — > [8] 2 = 1000 — >• 0.0001 = [^ 2 (8)] 2 — >• /3 4 

9 — >• [9] 2 = 1001 — > 0.1001 = [0 2 (9)] 2 — ► /3 + /3 4 
10 — >• [10] 2 = 1010 — > 0.0101 = [</> 2 (10)] 2 — > /3 2 +/3 4 

e così via. 

Osservazioni 2.4.6. 1. La 1, 1-successione di punti {£™i} n eN presenta una 
forte analogia con la costruzione della successione di van der Corput, anche 
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se ad alcuni numeri espressi in forma binaria non corrisponde alcuna potenza 
di /3. Questo non ci deve meravigliare, se si considera da una parte che solo 
le rappresentazioni binarie corrispondenti alle sole composizioni possibili tra 
ipQ e ipi sono quelle a cui corrispondono potenze di (3, dall'altra che nella 
costruzione di van der Corput al passo n vengono divisi tutti i 2 n intervalli, e 
dunque si raddoppiano i punti (mentre per le L.S'-successioni si dividono solo 
gli intervalli lunghi). 



2. Infine, attraverso questa procedura "alla van der Corput' 
grado di determinare i punti della successione {Cì'iJneN svincolandoci dagli 
insiemi A™ x . Precisamente, i vantaggi legati all'introduzione di queste due 
funzioni ip e ipi sono due. Il primo è che per determinare i punti associa- 
ti all'rz-esima partizione p xi non è più necessario riordinare i punti di Pi^ 1 
perchè essi si ottengono scrivendo tutte le n-uple ordinate di elementi dell'in- 
sieme {0, 1} (che sono 2 n ), cancellando le sequenze di n cifre che presentano 
due cifre consecutive pari a 1. 
Questo fatto si vede bene in Af x , Af x e A* x . 

L'altro vantaggio è che, mentre nella Definizione |2.3.1 e nell'Esempio 2.3.3 



ogni insieme A™ : è ottenuto attraverso funzioni dipendenti da n, in que- 
sto caso esistono per tutti gli n solo le funzioni ip e ipi (indipendenti da 
n), e ogni punto si ottiene attraverso composizioni successive di queste due 
funzioni. 



3. La 1, 1-successione di punti {^"xjnGN può essere ottenuta utilizzando 
la rappresentazione binaria dei numeri naturali, invertendo l'ordine delle sue 
cifre (ed ottenendo, cioè, la rappresentazione binaria della funzione radice 
inversa) ed eliminando le rappresentazioni che presentano due cifre conse- 
cutive pari ad 1. Inoltre, i coefficienti delle potenze di (3 sono esattamente 
quelli che compaiono nella rappresentazione binaria di 02, come nel caso della 
successione di van der Corput. Notiamo, inoltre, che c'è una corrisponden- 
za biunivoca tra i coefficienti delle potenze di /3 e la sua rappresentazione 
attraverso composizioni successive delle funzioni ip e Se consideriamo, 
ad esempio, Af A , vediamo che Vo, i>i, 4>oi V^ooi, V'ioi, ^oooi, V'iooi corrispon- 
dono rispettivamente a 0.00, 0.10, 0.01, 0.001, 0.101, 0.0001, 0.1001. E le 
cancellazioni delle composizioni non lecite (e cioè tutte quelle che usano V'i.i) 
corrispondono esattamente (nell'ordine giusto) alle cancellazioni dei valori di 
[0 2 (')]2 le cui sequenze di cifre contengono la coppia 11. 



Passiamo ora agli altri due esempi, prima di fornire la regola generale. 
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Esempio 2.4.7. Consideriamo la successione di punti {£2,i}neN partendo 
dalla successione di partizioni {p£ i}neN ottenuta per L — 2, S — 1. Definia- 
mo tre funzioni sui sottointervalli di [0, 1] nel modo seguente 

V>o:[0,l[— >[O,0[ (2.9) 
x i — y (3x , 

V>i:[0,l[— >[/9,2/9[ (2.10) 
xt-+ /3x + /3 

e 

^2 : [0,/3[ — ^ [2/5, 1[ (2.11) 
x H> f3x + 2/3 . 

In questo caso, ovviamente, le uniche composizioni proibite sono ip2 ° V'i e 
■02 ° V ; 2- Valutiamo ^o, V'i e ^2 nello zero e nuovamente otteniamo l'insieme 
dei punti associati alla prima partizione ply. 

Aii = (eii.^i.Si) = (^o(0),Vi(0),^(0)) = (0, /3, 2/3) . 

Calcolando, ove possibile, le funzioni (2.9), (2.10) e (2.11) nei punti ottenuti, 
ossia componendo tra loro queste funzioni, come nel caso della successione 
di punti di Kakutani-Fibonacci, si ottiene l'insieme dei punti relativi alla 
partizione successiva, cioè 

A 2 - (V £ 2 f 3 f 4 f 5 f 6 f \ 

iv 2,l ~~ l S2,l> S2,l> S2,l> S2,1>S2,1> S2,l> S2,l J 

= (^(O),^^),^^),^!^),^!^),^^),^^)) 

= (^oo(O), Vio(0), ^o(O), Voi(0), Vn(0), Vo 2 (0), Vi 2 (0)) 
= (0,/3,2/3,/3 2 ,/3 + /3 2 ,2/3 2 ,/3 + 2/3 2 ) . 

Anche in questo caso, come nell'esempio precedente, si ha ipio = ipi(0) per 
ogni « = 0,1,2. Iterando un'ultima volta questa procedura si ha 

A3 _ (ti £2 t3 £4 £5 t6 SI f 8 
1V 2,1 — l S2,l) S2,l; S2,l) S2,l; S2,l; S2,l) S2,l) S2,l; 
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t% tlO tll £12 tl5 d-16 tl7\ 
S2,l> S2,l> S2,l> S2,l) S2,l) S2,l) S2,l ) 

= (^0(0)^1(0)^2(0)^01(0)^11(0)^02(0), 

^12(0), Vooi(O), Vioi(O), Von(O), Vm(0), 
^002(0), ^102(0), ^012(0), ^112(0)) 

= (Vooo(O), Vioo(O), ^200(0), Voio(O), ^no(O), ^020(0), 
^120(0), Vooi(O), Vioi(O), ^201(0), V011 (0), (0), 
V002 (0) , V102 (0) , V202 (2) , V012 (0) , V112 (0)) 

= (o,/3,2/3,/3 2 ,/3 + /3 2 ,2/3 2 ,/3 + 2/3 2 ,/3 3 , 
f3 + /3 3 , 2/3 + /3 3 , /3 2 + /3 3 , /3 + /3 2 + /3 3 , 
2/3 3 , /3 + 2/3 3 , 2/3 + 2/3 3 , /3 2 + 2/3 3 , /3 + /3 2 + 2/3 3 ) , 

e così via si ottiene {C^ilneN- 

I punti ottenuti da composizioni ricorsive di queste funzioni risultano essere 
disposti secondo l'ordine giusto. 

Vogliamo vedere quale relazione intercorre tra queste funzioni e la rappresen- 
tazione dei numeri naturali in base 3 e se è lecito supporre che le composizioni 
non consentite corrispondano a rappresentazioni di numeri naturali che non 
danno luogo ad alcun punto della successione. Infatti 




1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



[0]a 
[l]a 
[2]a 
[3]a 
[4]a 
[5]a 
[6]a 
[7]a 
[8]a 
[9]a 
[10]a 



00 
01 
02 
10 
11 
12 
20 
21 
22 
100 
101 



0.00 
0.10 
0.20 
0.01 
0.11 
0.21 
0.02 
0.12 
0.22 
0.001 
0.101 



[M0)h 
[0 3 (i)] 3 

3 (2)] 3 

[03(3)] 3 

[0 3 (4)] 3 
[0 3 (5)] 3 
[0 3 (6)] 3 
[0 3 (7)] 3 
lM8)h 
[0 3 (9)] 3 
[0 3 (1O)] 3 





/? 
2/3 

/? 2 

/? + /3 2 
nessun /3 
2/3 2 

/? + 2/3 2 
nessun /3 



e cosi via. 



Osservazione 2.4.8. La successione di punti {£2 i}neN può essere ottenuta 
dalla rappresentazione ternaria dei numeri naturali, invertendo l'ordine delle 
sue cifre ed eliminando le rappresentazioni che presentano due cifre consecu- 
tive del tipo 21, 22. Le potenze di /3 possono essere calcolate come nel caso 
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della successione di van der Corput, tenendo conto dei coefficienti moltipli- 
cativi; precisamente, le cifre nella rappresentazione di 03 (n) in base 3 sono 
proprio i coefficienti di (3. 

Esempio 2.4.9. Consideriamo come ultimo esempio la successione di punti 
{£" 2 }neN e determiniamo la famiglia di funzioni per la successione di parti- 
zioni p™ 2 , ottenuta mediante successivi p-raffinamenti della partizione banale 
cu, con p = {[0, (3[, [(3, (3 + (3 2 [, [(3 + (3 2 , 1[}. Definiamo le seguenti funzioni: 

^ :[0,1[— >[0,/9[ (2.12) 
x \-> j3x , 

+ 2 [ (2.13) 

x !->■ f3x + (3 

e 

^ 2 : [0,/5[ — > [/3 + /5 2 ,l[ (2.14) 
x^ f3x + (3 + f3 2 . 

Anche in questo caso analizziamo le composizioni possibili tra queste funzioni. 
Osservato che ben quattro delle nove possibili composizioni non sono lecite, 
e cioè ipi o ipioi/j 2 , Ì J 2°Ì J i e Ì J 2°Ì J 2, passiamo subito alla determinazione 
dei primi tre insiemi di punti relativi alle partizioni p\ 2 , p\ 2 , p\ 2 . 

Al, 2 = = (Vo(0),Vi(0),V 2 (0)) = (0, 0, p + /3 2 ) . 

Componendo le funzioni in maniera lecita otteniamo 
A 2 - (V £ 2 f 3 f 4 f 5 ì 

Ji l,2 ~~ l Sl,2) Sl,2> S1,2)S1,2) Sl,2^ 

= (^ (0),Vi(0),V2(0),Voi(0),Vo2(0)) 
= (Voo(0), ^io(O), V 2 o(0), Voi(0), Vo 2 (0)) 
= ( ,/3,/3 + /3 2 ,/3 2 ,/3 2 + /3 3 ) . 

Iterando nuovamente questo procedimento ai punti ottenuti si ha 

a 3 _ /Vi £2 t3 ^4 t5 t6 t-7 d-8 t9 ^10 tll \ 
Yi l,2 — l S1,2)S1,2) Sl,2) S1,2)S1,2) Sl,2) S1,2)S1,2>S1,2) Sl,2) Sl,2 J 



2.4 Procedimento "alla van der Corput" 



40 



= (Vo (0) , Vi (0) , V>2 (0) , Voi (0) , V02 (0) , Vooi (0) , 

Vioi (0) , V201 (0) , V002 (0) , V102 (0) , V202 (0) ) 

= (Vooo (0) , V100 (0) , V200 (0) , V010 (0) , V020 (0) , V001 (0) , 

Vioi(0), (0), Voo 2 (0), $102(0), $202(0)) 
= (0,/3,/3 + /3 2 ,/3 2 ,/3 2 + /3 3 ,/3 3 ,/3 + /3 3 , 

/3 + /3 2 + /3 3 , /3 3 + (3\ (3 + /3 3 + /3 4 , /3 + /3 2 + /3 3 + /3 4 ) . 

Procedendo in questo modo si ottiene la successione di punti {Ci^}- 
Osserviamo che anche in questo caso si tratta di considerare la rappresenta- 
zione dei numeri naturali in base 3. Vedremo che ci sono due cancellazioni 
in più da effettuare rispetto al caso L — 2, S — 1, relative alle composizioni 
non possibili, ma l'osservazione cruciale riguarderà i coefficienti di (3. 
Vediamo, a tal scopo, cosa succede partendo dalla rappresentazione in base 3 
dei primi numeri naturali, così come abbiamo fatto nei due esempi precedenti. 

— > [0] 3 = 00 — > 0.00 = [fo(o)] 3 — > 

1— ► [l]s = 01— > 0.10 = [0 3 (1)] 3 — > p 

2^ [2] 3 = 02^ 0.20 = 3 (2)] 3 — > (1 + ^)13 

3^ [3] 3 = 10— ► 0.01 = [0 3 (3)] 3 — ► /? 2 

4 — ► [4] 3 = 11 — >■ 0.11 = [<fo(4)] 3 — >■ nessun /3 

5 — >• [5] 3 = 12 — > 0.21 = [<fo(5)] 3 — >■ nessun /3 
6^ [6] 3 = 20^ 0.02 = [0 3 (6)] 3 — ► (1 + /3)/3 2 
7^ [7] 3 = 21—». 0.12 = [0 3 (7)] 3 ->■ nessun /3 
8 — ► [8] 3 = 22 — > 0.22 = [<A 3 (8)] 3 — > nessun fi 
9^ [9] 3 = 100^ 0.001 = [0 3 (9)] 3 — ► /? 3 

10-». [10] 3 = 101— > 0.101 = [<fe(10)] 3 — > /? + /3 3 

e così via. 

Osservazione 2.4.10. Come nei casi precedenti, anche in quest'ultimo la 
successione di punti 2 }neN può essere ottenuta dalla rappresentazione in 
base 3 dei numeri naturali, invertendo l'ordine delle sue cifre ed eliminando le 
rappresentazioni corrispondenti a composizioni non lecite. Ciò che rende que- 
sto esempio particolarmente interessante è l'espressione dei coefficienti delle 
potenze di (3. Infatti, mentre per gli esempi corrispondenti alle successioni di 
punti {£™i}neN e {^2,i}neN, rispettivamente, questi erano proprio le cifre di 
02 (n) e 03 (n), rispettivamente, in base 2 e 3, rispettivamente, in questo caso 
ciò è falso per la cifra 2. Al posto di 2 bisogna sostituire l'espressione 1 + f3. 
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Siamo pronti per considerare il caso generale per L, S G N. 



Il caso generale s } ne N 

Supponiamo di voler suddividere l'intervallo unitario in L intervalli lunghi di 
ampiezza (5 ed S intervalli corti di ampiezza /3 2 , con Lf3 + S(3 2 = 1. Definiamo 
le seguenti funzioni 

^ : [0,1[ — ^ [j/3, (j + l)/3[ (2.15) 
x i — y f3x + j(3 

con j = 0, . . . , L — 1 e 

4>L,k ■ [0, /?[— > [£/3 + L/5 + (fc + l)/3 2 [ (2.16) 
x ^ /3x + L/3 + k/3 2 , 

con = 0, . . . , S — 1. 

Si vede immediatamente che gli insiemi di punti s sono ottenuti dalla 
composizione di queste funzioni, ma ovviamente non tutte le composizioni 
sono possibili. Precisamente, le composizioni non ammissibili sono quelle 
del tipo ipL^ ° i>j per ogni valore di k G {0, . . . , S — 1} e per ogni valore di 
j G {1, . . . , L - 1} e ìp LM o tp Ltka per ogni k u k s . 

Allora, valutando queste due famiglie di funzioni nello zero, troviamo l'insie- 
me 

Ai,s = {Vo(0), Vi(0), • • • , ^l-i(O), ^,o(0), ^,i(0) . . . , Vi,s-i(0)} 

dei primi L + S punti di ^jneN- Componendo in maniera lecita le funzio- 
ni, si ottengono tutti gli altri insiemi, come abbiamo già visto nei tre esempi 
precedenti nei casi particolari delle LS'-successioni di punti, con L = S = 1, 
L = 2,S = leL = l,S = 2. 

Anche nel caso generale, ovviamente, esiste uno stretto legame tra queste 
funzioni e le potenze di (3 (e la rappresentazione dei numeri in base L + S). 
Illustriamo qui di seguito la natura di tale legame. 

Innanzitutto, per rappresentare un numero in base L + S abbiamo bisogno 
delle cifre da a L + 5 - 1, ossia dei resti nella divisione di un generico 
numero naturale per L + S. Poi, dopo aver espresso ciascun numero naturale 
in base L + S, si inverte tale rappresentazione (ottenendo, così, la rappresen- 
tazione in base L + S della funzione radice inversa calcolata in tale numero) 
e si eliminano le sequenze di cifre che contengono al loro interno gruppi di ci- 



fre consecutive corrispondenti a composizioni non lecite delle funzioni (2.15) 



e (2.16). Come negli esempi precedenti, queste cifre rappresentano (nella 
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maniera appresso specificata) i "coefficienti" delle combinazioni lineari delle 
potenze di (3 che ci danno i punti della successione. Precisamente, le cifre 
corrispondenti a valori compresi tra ed L sono proprio i coefficienti di /3, 
mentre per quelle comprese tra L + ledL + S — 1 vale la seguente relazione: 



k^- L + k(3 



(2.17) 



con k — 1, . . . , S — 1, e cioè alla cifra L + k si sostituisce L + h/3. 
Chiaramente, per L — S — 1, L — 2eS — leL — leS — 2 otteniamo 
esattamente i risultati presentati gli esempi |2.4.5| |2.4.7| e |2.4.9 



Osservazione 2.4.11. Sia la Definizione 12.3.11 che il metodo che abbiamo 
appena esposto nel caso generale utilizzano ciascuna una propria famiglia 
di funzioni da applicare in maniera ricorsiva ai punti degli insiemi prece- 
dentemente ordinati A^g 1 per la determinazione degli insiemi di punti A£ s . 
Tuttavia, sebbene a prima vista le due definizioni possano sembrare sostan- 



zialmente equiparabili, il secondo metodo che usa le funzioni (2.15) e (2.16) 



risulta particolarmente efficiente per via di due considerazioni che sono la na- 
turale generalizzazione di quelle esposte nelle Osservazioni 2.4.6 La prima 



riguarda la natura delle funzioni (2.15) e (2.16), che in questo caso vengono 



definite all'inizio e quindi non più ridefinite per ogni singolo blocco, a diffe- 



renza delle funzioni (2.5), che dipendono da n. La seconda considerazione è 



forse più importante della prima e dipende dal legame stretto che intercor- 
re tra queste funzioni e la rappresentazione in base L + S, come abbiamo 
evidenziato qui sopra: infatti, gli insiemi di punti A£ $ vengono determinati 
indipendentemente dall'insieme A^ 1 in quanto questo calcolo viene tradotto 
in termini di rappresentazione dei numeri naturali rispetto alla base L + S. 
Questo fa sì che se vogliamo calcolare i primi N punti di {££ s}neNi con 
t n -i < N < t n , possiamo farlo senza dover necessariamente calcolare tutti i 



punti di A^ s , perchè non si tratta più di applicare le funzioni (|2.5|) all'in 



sieme A™g, ma di scrivere i primi N punti nella base L + S, ovviamente al 
netto delle opportune cancellazioni, così come spiegato in precedenza. 

Ora che abbiamo schematizzato il nostro ragionamento, possiamo passare 
alla presentazione dell'algoritmo per calcolare le successioni di punti {££ s }. 



2.5 Algoritmo 

Concludiamo questo capitolo con la descrizione dell'algoritmo veloce, discus- 
so nel paragrafo precedente. Esso sarà scritto utilizzando uno pseudocodice, 
ossia mediante un linguaggio di programmazione fittizio, non direttamente 
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compilabile, allo scopo di renderlo più leggibile e comprensibile. 

In seguito, per determinare le successioni di punti {££ s}neN, si è scelto di 

tradurre l'algoritmo in un linguaggio implementabile dal software Matlab. 

Procediamo alla descrizione dell'algoritmo in tre passi, la cui parte princi- 
pale sarà costituita dal corpo centrale. 

1. 

Algorithm: base Jnversa 
Input: (n G N, b = L + S) 
Output: [n]i 

Data: int j = [\og(n)/log(b)\ + 1 

int a[j] = 

begin 

for i — 1 — >■ j do 
a[i] <— mod(n, b) 
n <- Lf J 
end 
end 

Result: a 

Questa seconda sequenza di comandi, alquanto complessa visto il numero 
di cicli innestati, richiama anche l'algoritmo precedente, che abbiamo chia- 
mato funzione baseJnversa e rappresenta la prima parte vera e propria del 
nostro algoritmo, in cui vengono effettuate le cancellazioni dei vettori a che 
presentano cifre consecutive corrispondenti alle composizioni non ammissibi- 
li. Il risultato di questo secondo algoritmo sarà una matrice A, le cui righe 
sono i vettori a che restano dopo aver effettuato le cancellazioni (che cor- 
rispondono ai coefficienti delle potenze di (3 nella maniera che vedremo nel 
terzo algoritmo). 

E dunque sufficiente fornire in input solo i tre valori di n, L, S. 
Notiamo innanzitutto che, poiché non tutte le rappresentazioni dei numeri 
naturali in base L + S corrispondono a punti della successione, non sappiamo 
a priori a quale sequenza di cifre corrisponderà l'n-esima riga della matrice 
(che non può, evidentemente, corrispondere al numero naturale n). A causa 
di questa incertezza sul numero di righe, è dunque necessario inizializzare un 
vettore riga A. Per la stessa ragione abbiamo poi scelto che questo vettore 
riga abbia dimensione 100, anche se questo numero potrebbe essere ridotto, 
dal momento che con 100 cifre possiamo esprimere un numero estremamente 
grande in qualsiasi base. Inoltre, poiché il nostro scopo è quello di ottenere 
la matrice A, aggiungiamo passo passo altre righe al vettore A, le cui com- 
ponenti sono tutte nulle per rappresentare lo zero. Queste righe successive 
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saranno identificate con un vettore ausiliario, che indichiamo con c. Ciò che 
faremo sarà controllare ad ogni passo il numero di righe della matrice A, che 
denotiamo con r. Fintanto che il numero di righe è minore del numero fissato 
n, prenderemo in esame la rappresentazione del numero successivo, indicato 
con i. A questo punto bisogna controllare se questo numero ha una rappre- 
sentazione che contiene delle cifre corrispondenti a composizioni di funzioni 
ammissibili oppure no, e per farlo utilizziamo una variabile ausiliaria d, il cui 
valore sarà o 1 nel caso in cui il numero debba essere inserito nella matri- 
ce o no, rispettivamente. La lunghezza del vettore a è variabile e dipende 
dal numero che stiamo rappresentando. Ogni volta che rappresentiamo un 
numero naturale, otteniamo a con l'algoritmo 1, immagazziniamo le sue com- 
ponenti ed effettuiamo un controllo su esse. Se due componenti consecutive 
di a coincidono con le composizioni non lecite, poniamo d — d + 1, altrimenti 
non si fa nulla. Nel primo caso, incrementiamo semplicemente il valore di 
i, mentre nel secondo, in cui il valore di d è ancora impostato sullo zero, 
memorizziamo nelle prime cifre di c (numero pari alla lunghezza di a) gli 
stessi valori contenuti in a e controlliamo se sono presenti valori maggiori di 
L. In tal caso riproponiamo il ragionamento fatto nel Paragrafo 2.17 nel cui 
caso sostituiamo a ciascuna cifra maggiore di L l'espressione (2.17). Infine 
incrementiamo il valore dell'intero positivo successivo % e aggiungiamo la riga 
c alla matrice A. Ora, il numero di righe r della matrice viene calcolato nuo- 
vamente e si ripetono così tutti i cicli. Alla fine, quando il numero di righe 
della matrice è n, si calcola il numero di cifre e necessario per rappresentare 
l'ultimo elemento i, in modo che la matrice finale A sia ripulita dalle colonne 
superflue, ossia A sarà la matrice costituita dalle prime r righe ed e colonne 
della matrice di partenza. 

3. 

Algorithm: Algoritmo 
Input: (L, S,neN) 
Output: {C L ,s} 

Data: A = cancellazioni(L, S, n) 

a _ -L+VL 2 +4S 

begin 

[•, c] = size(A) 
for k = 1 — > c do 
| s(k) = f3 k 
end 

B = As' 
end 

Result: B 
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Infine quest'ultimo algoritmo restituisce esattamente il numero di punti 
da noi indicato, ossia, fissato n G N, esso restituirà proprio i primi n termini 
della successione di punti {££ 5} n eN- Esso prende in input la matrice A, cal- 
cola il numero delle sue colonne e lo indica con c. Poi costruisce un vettore 
s di c componenti, ciascuna delle quali è una potenza crescente di (3. Quindi 
il nostro insieme di punti s> ££ 5, • • • , £™ 5) sarà il prodotto della matrice 
A con il vettore s. 



Poiché, come abbiamo già detto, il precedente algoritmo è stato imple- 
mentato con Matlab (che è un software numerico) la risposta che otteniamo è 
il valore numerico delle potenze di (5 corrispondenti. E ovviamente possibile 
anche una risposta simbolica da parte del programma, ma abbiamo scelto 
la prima per due ragioni: una puramente tecnica riguarda i grafici, che non 
sarebbero stati possibili attraverso un parametro simbolico, l'altra riguarda 
la velocità di risposta del computer, che impiega molto più tempo per elabo- 
rare una risposta simbolica. 



Riportiamo le risposte del programma nei tre Esempi 2.3.3 considerati in 
precedenza. 

1. I primi 8 punti della successione {^"ijneN; corrispondenti a Af^ sono 

0, 0.6180, 0.3820, 0.2361, 0.8541, 0.1459, 0.7639, 0.5279 . 

2. I primi 17 punti della successione {£2 i}neN) corrispondenti a Af t sono 

0, 0.4142, 0.8284, 0.1716, 0.5858, 0.3431, 0.7574, 0.0711, 0.4853, 
0.8995, 0.2426, 0.6569, 0.1421, 0.5563, 0.9706, 0.3137, 0.7279 . 

3. Infine, i primi 11 punti della successione {£i i2 }neN, corrispondenti a Af 2 
sono 

0, 0.5000, 0.7500, 0.2500, 0.3750, 0.1250, 
0.6250, 0.8750, 0.1875, 0.6875, 0.9375 . 
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Capitolo 3 

LS'-successioni di punti in 

[0,l[x[0,l[ 

In questo terzo ed ultimo capitolo estenderemo la definizione di LS'-successioni 
di punti e l'algoritmo presentato nel capitolo precedente al quadrato di lato 
unitario I 2 = [0, l[x[0, 1[. Quest'ultimo ci sarà fondamentale per la visualiz- 
zazione, mediante dei grafici, di alcuni dei risultati ottenuti. 
Per far ciò, prenderemo in esame le successioni introdotte da van der Corput 
in I 2 [11] e da Halton in P |5j, restringendo la nostra attenzione al caso 
bidimensionale. 



3.1 Successioni di van der Corput e di Halton 

Richiamiamo preliminarmente alcuni concetti esposti nei capitoli precedenti. 
Ricordiamo che se b è un intero positivo, allora per ogni numero naturale 
n = a + ai 6 + . . . + aM-i& M_1 + a>Aib M l'espressione 

[njj, = dM a M—l ■ ■ • &2 a l a 

indica la rappresentazione di n in base b, dove M = |_log b nJ è il numero di 
cifre necessarie per esprimerla. 

Come abbiamo già avuto modo di osservare nella Definizione 1.4.9 invertendo 
l'ordine delle cifre a^, con % e {0, . . . , M}, si ottiene il numero compreso tra 
e 1 detto radice inversa di n, e cioè 

(f) b (n) = aob- 1 + a x b- 2 + ... + a M b~ M ~ l , (3.1) 

la cui rappresentazione in base b è [0&(n)]6 = 0.a ai . . . a,M-\ a M- 
Siamo adesso pronti per ricordare la definizione di alcune importanti succes- 
sioni che fanno uso di tale notazione e che in qualche modo generalizzano in 
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spazi multidimensionali la successione di van der Corput già presentata nella 
Definizione 11.4.111 

La prima interessante generalizzazione in dimensione 2 risale allo stesso van 
der Corput che in [11] ha introdotto la seguente successione finita. 

Definizione 3.1.1 (Successione di van der Corput). La successione di van 
der Corput di ordine iVeN nel quadrato I 2 è definita come 

x n = (-^r,0 2 (n)J con n = l,...,N. (3.2) 

L'immagine successiva riporta la rappresentazione grafica della successio- 
ne (5000' M n )) per n = 1, . . . , 5000. 




Figura 3.1: Successione di van der Corput ( , <fe(^)) per 1 < n < 5000 



In seguito, fu il matematico inglese J. M.Hammersley che in [Bj propose una 
versione s-dimensionale della definizione precedente. Egli, infatti, prende i 
primi s — 1 numeri primi, che indica con b\, 62, • • • , &s-i- Considera poi la rap- 
presentazione dei primi N numeri naturali rispetto a queste basi e il processo 
di riflessione di tali espressioni in modo da ottenere frazioni nell'intervallo 
[0,1] del tipo (|3l]). 
Si ha dunque la seguente 
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Definizione 3.1.2 (Successione di Hammersley) . La successione di Ham- 
mersley di ordine N e N nell'ipercubo P è la successione s-dimensionale 
definita come 

*n= (S>&i( T 0>&a( T 0>--->&.- 1 (w)) con n = 1,...,JV , (3.3) 
dove 6 S _! sono i primi s — 1 numeri primi e le 0&. sono le funzioni 



radice inversa che compaiono nella Definizione |1.4.9 



Osserviamo che la successione di Hammersley è un insieme di iV pun- 
ti nell'ipercubo I s e che quindi tale insieme non può essere esteso ad una 
successione infinita: infatti, fissato il numero N, non si possono aggiungere 
ulteriori punti a questa sequenza, a meno di non aumentare il valore di N, 
ma ciò comporta che tutti i termini vengano ricalcolati di conseguenza. 

La successione di van der Corput fu generalizzata in un'altra direzione da 
Halton nel modo seguente. 

Definizione 3.1.3 (Successione di Halton). Definiamo successione di Halton 
s-dimensionale (x) n6 N la successione definita sull'ipercubo I s come 

x n = (0 6l (n),0 62 (7i),...,0 6s (n)) per ogni n E N , (3.4) 

dove b%, & 2 , ■ ■ ■ , b s sono interi positivi a due a due coprimi. 

È immediato osservare che per la successione di Halton la dipendenza 



da N è stata eliminata. Nella Figura |3.2| riportiamo il grafico dei primi 
5000 punti della successione di Halton bidimensionale, relativi ai due numeri 
coprimi b\ = 2 e b 2 = 3. 

Ricordiamo che le successioni di punti in I s per cui la discrepanza risulta 
essere uguale a C(log s N/N) per iV — > oo sono dette successioni a bassa 
discrepanza. 



Lo stesso Halton in [S] ha dimostrato che la (3.4) è una successione a bassa 
discrepanza, cioè vale il seguente risultato 

Proposizione 3.1.4. La discrepanza della successione di Halton s-dimensio- 
nale (x n ) nf z N soddisfa 

D N (xi, x N ) < — 

per ogni n e N. 

Si può infatti osservare che i punti della successione { (y9 2 (n), ^(tì)) } 



nei 



illustrati nella Figura 3.2 , si distribuiscono in modo da ricoprire il quadrato in 



maniera uniforme, come ci si aspetta da una successione a bassa discrepanza. 
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0.2 0.4 0.6 0.8 1 



Figura 3.2: Successione di Halton {((^(n), 03(>i))}neN per 1 < n < 5000 

Anche le successioni di van der Corput e Hammersley sono esempi di suc- 
cessioni a bassa discrepanza, ed entrambi gli autori presentarono risultati 
importanti riguardanti la stima della discrepanza. 

Le successioni a bassa discrepanza rivestono un ruolo molto importante nelle 
applicazioni perchè si tratta di successioni utilizzate nei metodi Quasi- Monte 
Carlo per il calcolo approssimato di integrali. 

Sappiamo che esistono vari metodi di integrazione numerica per il calcolo 
di integrali di funzioni di cui non riusciamo a determinare analiticamente la 
primitiva. Questi metodi approssimano il valore dell'integrale della funzione 
/ con una combinazione lineare di valori assunti dalla /, ad esempio, nel caso 
unidimensionale 

N 

/ f{x)dx « y^jfjXj) . 
J 1 n=l 

Ovviamente differenti formule di quadratura corrispondono a differenti scel- 
te dei punti Xi e dei pesi Wi. Il metodo di integrazione Monte Carlo si basa 
sull'utilizzo di numeri Xi casuali, con pesi Wi = per i — 1, . . . , N. 
Se, ad esempio, vogliamo determinare l'integrale di una funzione / sul qua- 
drato I 2 , consideriamo i vettori casuali xi, X2, . . . , x^v e approssimiamo il 
valore dell'integrale con la media dei valori assunti dalla funzione in detti 
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punti, ossia 



/ /(x)dx 



1 



N 



Jl 2 



N 



n=l 



La caratteristica peculiare di tale metodo consiste nel fatto che l'errore asso- 
luto 



ossia l'errore commesso nella stima dell'integrale attraverso la media dei valo- 
ri della / negli N punti Xi, decresce in media come in virtù del Teorema 
del Limite Centrale, indipendentemente dalla dimensione dello spazio. 
Il metodo Monte Carlo può dunque sembrare alquanto insoddisfacente se con- 
frontato con altri metodi di quadratura noti per l'approssimazione di integrali 
di funzioni in una variabile come il metodo dei trapezi, di Cavalieri-Simpson 
o di Gauss, che convergono molto più rapidamente, nel senso che l'errore 
decresce più velocemente. 

Quando però si passa alla stima di integrali multidimensionali la situazione 
cambia notevolmente, perchè mentre l'errore con la stima fornita dal metodo 
Monte Carlo risulta sempre dello stesso ordine, con gli altri metodi l'ordine 
di convergenza decresce al crescere della dimensione dello spazio. 
Poiché la generazione di numeri o vettori casuali non è agevole, specialmente 
se N è grande, a partire da von Neumann, e quindi dai primi calcolatori, è 
stata sviluppata l'idea di considerare successioni deterministiche che simulino 
la casualità. Nei vari linguaggi di programmazione c'è di solito un comando 
denominato rand, o simili, che produce successioni di numeri di [0, 1[ che 
vari test statistici non distinguono da successioni casuali. Tali successioni si 
chiamano pseudo-casuali. 

Successivamente è nata l'idea di rinunciare ad imitare la casualità delle suc- 
cessioni e di considerare invece delle successioni in un certo senso ottimali 
per l'uso che si ha in mente, cioè per il calcolo degli integrali. Tali successioni 
si chiamano quasi-casuali (anche se il nome può ingannare, perchè non han- 
no nulla di casuale), che poi sono le nostre successioni a bassa discrepanza. 
Queste successioni danno luogo al metodo Quasi-Monte Carlo. 
Questa scelta pone rimedio a due aspetti: il primo riguarda la difficoltà di 
scegliere dei punti in maniera del tutto casuale, il secondo invece riguarda 
l'accuratezza nella convergenza dei metodi. Infatti sappiamo che una succes- 
sione di punti in / è a bassa discrepanza se essa è asintoticamente equivalente 
a log N/N e quindi l'utilizzo di questo secondo metodo garantisce che la stima 
dell'integrale convergerà al suo valore reale come ^jP-, ossia più velocemente 




di 



ì 
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Nei due articoli di Hammersley e Halton viene presa in esame la stima 
dell'integrale di una funzione / continua su I s , e cioè di 



dxi 



dx s f(xx, ...,x s ) 



(3.5) 



Hammersley in [6j afferma che, sebbene nei metodi numerici per il calco- 
lo approssimato di integrali si utilizzi, come abbiamo già visto, la formula 
Y2i w if( x i)i i n assenza di altre informazioni sulla /, non possiamo far altro 
che considerare i pesi Wi tutti uguali. Poiché per ipotesi la / è necessariamen- 
te integrabile, essa può essere approssimata da una successione di funzioni 
semplici integrabili. Inoltre, poiché gli insiemi su cui la / è costante si pos- 
sono costruire mediante unione ed intersezione di intervalli, è ragionevole 



testare l'efficienza del metodo di approssimazione dell'integrale (3.5) suppo- 
nendo che la / sia la funzione caratteristica di un sottointervallo J* = [0, b[ 
contenuto nell'ipercubo unitario I s . 



Van der Corput, nel costruire la successione finita (3.2) sul quadrato, ha inol- 
tre dimostrato che esiste una costante C tale che, posto x n = (-^, (pz{n)) per 
ogni 1 < n < N, si ha 



sup 

6i,6 2 e[0,l[ 



A' 



N{b x b 2 



n=l 



<C\nN 



dove b = ibi, b 2 ). 

Hammersley, dopo aver introdotto la successione finita (3.3), si chiede in [5] 
se è vero che 



sup 

j*ci s 



N 



J2xj*M - N\(J* 



n=l 



sia maggiorato da un multiplo di una qualche potenza di logiV. 

Halton dimostra in [5] che per la successione di Hammersley (x n )i< n <jv 



definita da (3.3) vale la seguente stima: 



sup 

j*ci» 



N 



^ Xj ,(x n )-iVÀ(J* 



n=l 



<C s (lnN) 



s-l 



per ogni neN, con C s esplicita costante positiva. 



3.2 LS'-successioni "alla van der Corput" 

In questo capitolo presentiamo una prima estensione al caso bidimensiona- 
le delle LS'-successioni di punti. Vogliamo presentare una loro variante in 
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maniera simile a quanto fatto da van der Corput nel passaggio dalla succes- 
sione unidimensionale {02(^)}neN alla successione finita 4>2(n)) 1<n<N nel 
quadrato. 

Definizione 3.2.1. Una LS -successione "alla van der Corput" di ordine 
ÌVgN nel quadrato I 2 è definita come 

n \ 

-,Cl,s) <xmn = l,...,N, (3.6) 

dove {££ s}n6N è una LS-successione. 

Ricordiamo brevemente la definizione di discrepanza per una successione 
in I 2 , già presentata nel più generale contesto multidimensionale in Par. 1.3, 
perchè in questo paragrafo ci occuperemo esclusivamente del caso bidimen- 
sionale. 

Definizione 3.2.2. Sia x 1; ...,xjv una successione finita nel quadrato I 2 . 
La discrepanza D N di questi punti è definita come 



D N = Djv(x!, . . . , xat) = sup 



1 N 



n=l 



dove J = [ai, &i[x [a 2 , b 2 [c I 2 varia in tutti i modi possibili, con < aj < 
h <l,j = l,2. 



Di seguito riportiamo alcuni grafici relativi alle successioni di tipo (3.6) 
che ci permettano di visualizzare la distribuzione delle coppie nel quadrato e 
di fare qualche ipotesi sulla loro distribuzione. In particolare, sarebbe inte- 
ressante scoprire se si tratta ancora di successioni uniformemente distribuite, 
così come accade con la successione di van der Corput nel passaggio dall'in- 
tervallo unitario al quadrato. 



Le immagini seguenti mostrano una distribuzione delle coppie che riveste il 
quadrato in maniera piuttosto uniforme. Non ci sono infatti zone in cui i 
punti si addensano nè zone di rarefazione. La distribuzione dei punti in que- 
sti ultimi due casi, inoltre, appare addirittura migliore di quella di Halton in 



Figura 3.2 e confrontabile con quella di van der Corput in Figura 3.1 
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Figura 3.4: Successione (5550, £21) per ^ — n — ^00 
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Figura 3.5: Successione (5550,^31) P er 1 < n < 5000 

Nel prossimo ed ultimo paragrafo introdurremo delle L5"-successioni infinite, 
prenderemo in esame alcuni significativi esempi e cercheremo di determinare 
quale successione ha il comportamento migliore. 

3.3 L^S-successioni "alla Halton" 

Volgiamo ora la nostra attenzione ad una seconda, ultima e più complessa 
generalizzazione delle L^-successioni di punti nel quadrato unitario. Esse 
vengono definite in questo contesto in modo naturale, in maniera analoga 
a quanto visto per le successioni di Halton. Infatti, l'algoritmo presentato 
nel capitolo precedente in Par. 2.5 ci permette di calcolare i punti della 
successione {££ s} n eN senza doverci preoccupare della posizione di tali punti 
in qualcuno dei blocchi A£ a differenza di quanto succedeva nella prima 
definizione di LS'-successioni di punti introdotte in [2]. 

Allora appare del tutto naturale considerare coppie di punti appartenenti a 
due LS- successioni (di punti). 

Definizione 3.3.1. Definiamo LS -successione "alla Halton" nel quadrato 
I 2 la successione {(££ l5l , ££ 2 <j 2 )}neN, dove {6£ l5l }„ 6N e {££ 2& }neN sono due 
LS'-successioni di punti. 
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A differenza della Definizione 3.2.1 con la quale abbiamo introdotto una 
successione finita di coppie di punti in I 2 , che varia al variare di N e N, in 
questo caso abbiamo definito una successione infinita. 

Non solo. Questa generalizzazione fa intravvedere due possibili generaliz- 
zazioni alle dimensioni s > 2. Potremmo infatti considerare s-uple di LS- 
successioni, oppure s-uple la cui prima componente sia la successione e 
le restanti s — 1 siano costituite da /^-successioni. 

Possiamo implementare l'algoritmo 2.5 separatamente per le successioni in 
ascissa e per quelle in ordinata e considerare poi le coppie di punti di tipo 
LS ottenute. In particolare abbiamo preso in esame, sia in ascissa che in 
ordinata, successioni di punti {^sì a Dassa discrepanza, ossia quelle relative 
alla relazione L/3 + S(3 2 = 1, con L > S. 

Qui di seguito riportiamo i grafici più interessanti tra quelli ottenuti in se- 
guito a queste implementazioni. 
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0.2 0.4 0,6 0.8 1 



Figura 3.7: Successione {(£3 ij£s i)}neN P er 1 < n < 5000 




0.2 0.4 0.6 0.8 1 



Figura 3.8: Successione {(^i, £J} i)}neN per 1 < n < 5000 



3.3 LS-successioni "alla Halton 



57 




Figura 3.10: Successione {(£31, C4i)}neN per 1 < n < 5000 
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Come si può osservare, si tratta di risultati molto diversi tra loro, seb- 
bene tutti discendano da successioni a bassa discrepanza su ciascun lato del 
quadrato. In particolare, risultano emblematici i primi due grafici, relativi 
alle successioni di coppie {(^i, £4 i)}neN e {(£3 1, £5 i)}neN, P er cui non s °l° 
non si conserva la bassa discrepanza, ma addirittura si è portati a ritenere 
che non si tratti di successioni uniformemente distribuite. 
Per quanto riguarda il primo grafico, la risposta a questa irregolarità potreb- 
be trovarsi nell'espressione dei (3 che risolvono le equazioni corrispondenti, e 
cioè: 

A/3 + /3 2 = 1 /3 = v / 5-2. 

I due valori risultano essere razionalmente dipendenti, ossia una loro combi- 
nazione lineare a coefficienti razionali dà come risultato un intero. Infatti 

- (V5-2) = 1 . 



È dunque lecito supporre che sussista un legame tra le due successioni di 
punti {£ii}neN e {£4i}«eN e fra le corrispondenti successioni di partizioni 
{Pl,lineN e {p^xIneN- 

Effettivamente, in questo caso la risposta è affermativa. 
Precisamente, vale il seguente risultato. 

Proposizione 3.3.2. Se t n è il numero totale degli intervalli di x e t' n il 
numero totale degli intervalli della partizione p\ x , allora si ha che 

t n = ^3n V er °9 n i weN . (3.7) 




Dimostrazione. Dimostriamo tale relazione per induzione. La (3.7) risulta 
vera per n = 1 poiché vale l'uguaglianza t x = 5 = t%. Supponiamo ora che la 
(3.7) sia vera per n — 1, ossia t n _ 1 = <3( n _i) e dimostriamo che è vera per n. 



alle differenze (2.4), e cioè 



Per quanto visto nel capitolo precedente, sappiamo che valgono le equazioni 

tn — t n -\ + t n _2 (3-8) 

4 = <-i + C-2 • (3-9) 
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Allora, utilizzando la prima di queste relazioni, si ha 

Hn = tZn-1 + ^3n-2 = Hn-2 + Hn~Z + Hn-2 = 
= ^3n-3 + ^3n-4 + ^3n-3 + ^3n-3 + ^3ra-4 = 
= 2(t 3rl _ 3 + £3^-4) + Ì3n-3 = 3t3„_3 + 2t 3n _ 5 + 2t 3n _ 6 = 
= 3t3„_3 + tzn-h + ^3n-3 — ^3n-4 + ^Hn-G = 
= 4t 3n _ 3 + Ì3„_6 , 

ma per ipotesi induttiva 4t 3 ( n _!) + £3^-2) = 4£^_ 1 + ossia la tesi. □ 

Attraverso la proposizione precedente abbiamo dimostrato che la succes- 
sione {t n } rìg N degli intervalli di {p2\} è una sottosuccessione della successione 



di Fibonacci {t n }neN richiamata in (3.8) 



A questo punto, vogliamo vedere se esiste una simile relazione anche tra {p^i} 
e {P5 x} e, più in generale, se è vero che, ad esempio, come nella Figura 3.7 



nella Figura |3.8[ la pessima distribuzione delle prime 5000 coppie della suc- 
cessione {(££,£5 £2',s')}neN nel quadrato dipenda dal fatto che la successione 
{^n}neN associata a {p r [ 5 } n eN sia una sottosuccessione di {t n } ne ^ associata 
a {p n L i s /}neN- 

Per determinare il valore che k deve assumere affinchè sia vera l'equazione 
t n = tkn per ogni n e N, osserviamo innanzitutto che essa deve essere vera 
anche per n — 1. Dunque, per risolvere il problema è sufficiente considerare 
l'equazione ottenuta per n — 1, e cioè 

t 1 = tk . 

Infatti, poiché = L' + S', il problema si riduce a trovare il k corrispondente 
a t k = V + S'. 

Nel caso delle successioni di partizioni {Pi^jneN e {p^ijneN bisogna osservare 
che £]_ = 4+1 = 5 = t*. che corrisponde a k = 3. 

Per quanto riguarda la coppia di successioni di partizioni {p^} e {p^} 
non esiste una relazione simile a quella trovata per {p^} e {p2i} perchè 
*! = V + S' = 5 + 1 = 6 non corrisponde ad alcun valore di £&. Dobbiamo 
dunque dedurre che non esiste alcun k che verifichi l'equazione t x = tk- 



Vediamo infine un confronto tra la successione di van der Corput di or- 
dine N, , 4>2( n )) j la successione di Halton {(0 2 (ra), 4>3(n))} nf z^, la LS- 
successione "alla van der Corput" di ordine N, (^,£31) e la LS'-successione 
"alla Halton" {(^3 x, ^JjneN presentate in questo capitolo, dal punto di vista 
dell'uniforme distribuzione. Per effettuare tale confronto prendiamo in esame 
la stima di un integrale doppio di una funzione continua di cui conosciamo il 
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valore esatto e vediamo quale successione fornisce la stima migliore. 
Sia f(x,y) = 2x + 3y 2 con (x,y) E I 2 . Allora naturalmente 



/ f(x, y)dxdy = 2x + 3y 2 dxdy = 2 . 
Jp Jp 

Nella tabella seguente riportiamo i valori di X^Li f( x i v)i 
con N = 100,500, 1000,2000. 



N 




{(0 2 (n),0 3 (n))} 


( Ni £3,1) 




Pseudo 
Random 


100 


1.9560 


1.9346 


1.9440 


1.958 


2.0592 


500 


1.9906 


1.9837 


1.9898 


1.9911 


1.9495 


1000 


1.9953 


1.9925 


1.9935 


1.9949 


2.0506 


2000 


1.9977 


1.9959 


1.9979 


1.9976 


1.9925 



Tabella 3.1: Stime dell'integrale di / 



Da questi calcoli risulta evidente che l'insieme di punti che meglio appros- 
sima il valore dell'integrale di / è quello corrispondente alla LS'-successione 
"alla van der Corput" di ordine N, (^,£3,1)- In particolare, questo risultato 
mostra un comportamento migliore di quest'ultima successione anche rispet- 
to alla successione di Halton. 

Per contro, si può osservare come sia poco soddisfacente l'approssimazione 
ottenuta con i numeri pseudo-casuali che presenta oscillazioni notevoli nel- 
l'errore di approssimazione. 

Visti i risultati sperimentali, riteniamo di poter formulare la seguente 

Congettura 3.3.3. Le LS -successioni "alla van der Corput" nel quadrato 
sono u.d.. 

Questo rimane comunque un problema aperto, come resta aperto il pro- 
blema della stima della discrepanza di queste successioni. 

Si tratta prevedibilmente di un problema difficile, poiché le analoghe dimo- 
strazioni per le successioni associate ai numeri primi (di van der Corput, 
Hammersley e Halton) hanno rappresentato e rappresentano delle pietre mi- 
liari nel campo della distribuzione uniforme. 

La congettura è d'altra parte supportata non solo dall'effetto visivo che ci 
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danno le figure di questo paragrafo, ma anche dai risultati numerici che ab- 
biamo presentato nella Tabella 3J. e che mostrano come i metodi svilup- 
pati in questo capitolo siano senza dubbio competitivi con quelli noti nella 
letteratura. 
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Capitolo 4 

Osservazioni conclusive 



La parte finale di questa tesi è stata dedicata alla costruzione di nuove suc- 
cessioni di punti nel quadrato [0, l[x [0, 1[ utilizzando uno strumento nuovo, 
introdotto recentemente in [2], e cioè le LS-successioni. Abbiamo, infatti, 
introdotto due classi numerabili di successioni di punti, che abbiamo chia- 
mato LS'-successioni "alla van der Corput" e LS- successioni "alla Halton". 
Successioni che appartengono alla prima classe sono di ordine N per ogni in- 
tero positivo N fissato (e, cioè, sono costituite da N elementi), mentre quelle 
appartenenti alla seconda classe sono infinite. 

Relativamente a questo problema abbiamo fatto molta sperimentazione 
sia grafica che numerica trovando dei risultati variegati e comunque interes- 
santi. Alcuni esperimenti hanno dato degli esiti francamente sorprendenti 
e che avranno bisogno di ulteriori approfondimenti teorici per spiegare certi 
fenomeni. 



Sono molto confortanti i risultati ottenuti considerando su una delle due 
coordinate la successione equidistribuita e sull'altra una LS-successione "alla 
van der Corput" (Figure 3.3, 3.4 e 3.5). Sia i risultati grafici che i risultati 



numerici ci hanno consentito di formulare la Congettura 3.3.3 che si può ri- 
tenere molto plausibile. 



Gli esperimenti che abbiamo fatto invece con coppie di LS'-successioni 
"alla Halton" hanno dato dei risultati contrastanti. Da un lato abbiamo ot- 
tenuto delle successioni nel quadrato che appaiono uniformemente distribuite 
sia con un numero basso di punti sia con un numero elevato di punti (Figure 



3.9 e 3.10) e sembrano competitive con le successioni di Halton per il calcolo 



dell'integrale di funzioni di due variabili. Ma abbiamo anche ottenuto delle 
successioni che palesemente non sono uniformemente distribuite e che danno 
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luogo a concentrazioni di punti in certe zone, accompagnate da altre zone 
prive di punti (Figure 3.6, 3.7 e 3.8). Un fenomeno analogo si verifica con le 
successioni di Halton associate a numeri che non siano coprimi. Ma mentre 
per le successioni di Halton la spiegazione era stata data da Halton stesso, 
per le L^-successioni manca ancora una spiegazione teorica per un comporta- 
mento tanto variegato. Una spiegazione (Proposizione 3.3.2) è stata trovata 
per l'accoppiamento {Ci i}neN e {^4i}neN, quello forse più sorprendente (Fi- 
gura 3.6). Per esso si verifica lo stesso fenomeno che si può osservare con 
l'accoppiamento di successioni di Halton quando si prendono basi non co- 
prime, e cioè una "risonanza' in corrispondenza al minimo comune multiplo 
delle due basi. Si sono però scoperti successivamente altri accoppiamenti 
che non sembrano dar luogo a successioni u.d. nel quadrato, per i quali tale 
risonanza non c'è. 



Questo fenomeno merita uno studio successivo, probabilmente non sem- 
plice, se si pensa alla lentezza del progresso della conoscenza in questo set- 
tore: ai lavori di van der Corput sono seguiti oltre vent'anni dopo quelli di 
Hammersley (che si era limitato a fare una congettura e degli esperimenti 
numerici) e di Halton, che ha fornito una base teorica per spiegare il buon 
comportamento delle successioni introdotte da Hammersley e ad introdurre 
a sua volta un'altra generalizzazione del procedimento di van der Corput. 

Un ulteriore filone di ricerca, qualora la congettura venga dimostrata, 
sarebbe indirizzato allo studio della discrepanza delle LS'-successioni u.d.. 
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